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ABSTRACT 

We study the ability of future CMB anisotropy experiments and redshift surveys to constrain a thirteen- 
dimensional parameterization of the adiabatic cold dark matter model. Each alone is unable to determine 
all parameters to high accuracy. However, considered together, one data set resolves the difficulties of the 
other, allowing certain degenerate parameters to be determined with far greater precision. We treat in 
detail the degeneracies involving the classical cosmological parameters, massive neutrinos, tensor-scalar 
ratio, bias, and reionization optical depth as well as how redshift surveys can resolve them. We discuss 
the opportunities for internal and external consistency checks on these measurements. Previous papers 
on parameter estimation have generally treated smaller parameter spaces; in direct comparisons to these 
works, we tend to find weaker constraints and suggest numerical explanations for the discrepancies. 

Subject headings: cosmology: theory - dark matter - large-scale structure of the universe - cosmic 
microwave background 



1. INTRODUCTION 



Current cosmological data have led astrophysicists to ex- 
plore a structure formation paradigm in which cold dark 
matter driven by adiabatic fluctuations lead s to the for- 
mation of g alaxies and CMB anisotropi cs (see Blumcntha] 



et al 



1984; Dodclson, Gates & Turner 1996 for reviews) 
Experiments planned for the next decade will be able to 
test this paradigm more stringently by searching for dis- 
tinctive features in t he power spectra o f CMB anisot ropies 
and polarization (see Hu et al. 1997 and Hu & White 1997b 
for reviews). If this framework is confirmed, then upcom- 
ing measurements, notably from the MAP[j and Planck^ 
satellites, will enable precision measurements of cosmolog- 
ical parameters such as the baryon fraction and matter- 



radiat ion ratio (Jungman et al. 1996a, b; Zaldarriaga et al 



1997; |Bond et al.| 1997; popcland et al.| 1998; ptoinpor k 
Efstathiou| 1998) 



As discussed by a number of authors, the details of the 
CMB power spectra contain a considerable amount of cos- 
mological information. However, this leverage is not com- 
plete; altering the model parameters in particular com- 
binations can yield power spectra that ar e observatio nally 
indist inguishable from a reference model (Bond et al. 1994, 



1997; |Zaldarriaga et al.| 1997). The presence of these so 
called degenerate directions means that a CMB data set 
will restrict the allowed models to a curve or surface in 
parameter space rather a point. 

The galaxy power spectrum of large redshift surveys 
such as the 2dF surveyf] and the Sloan Digital Sky Sur- 
vey (SDSS)^] provide a different window on the cosmologi- 
cal parameter space. Tak en alone, the resu lts are again 



plagued by degeneracies (Tegmark 199 7a; | Goldbergfc 



Strauss 1998; Hu et al. 1998, hereafter HET). However, 
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when combined with CMB data, the constraints can be 
highly complementary in that the degenerate directions of 
one lie along well-constrained directions for the other. A 
striking example involves the Hubble constant. By com- 
bining the data sets, features in the power spectrum can 
be measured in both real and redshift space, allowing the 
Hubble constant to be identified even though neither data 
set alone provid es a good constraint (Eisenstein et al. 1998, 
hereafter [EHr]). 

In this paper, we explore the details of how these data 
sets complement each other. In particular, we consider the 
constraints on cosmological parameters attainable by the 
statistical errors of this next generation of CMB experi- 
ments and redshift surveys. We identify physical mech- 
anisms by which the data sets resolve degeneracies and 
explore them through progressions of cosmological models 
in the baryon fraction, neutrino mass, and tensor contri- 
bution. The combination of CMB data and large-scale 



struct ure has been studie d with curr ent d ata (Scott et al 
1995; [Bond fc Jaffel 1 996; |Lineweaver| 1998; ICawiser & Silk 
1998; IWebster et al 



1998) as well as with future data in 
a smaller space of cosmolo gical paramet ers but with more 
general initial conditions (Wang et al. 1998). We study 
how reducing the cosmological parameter space affects the 
degeneracies. 

In addition to incorporating large-scale structure data, 
our treatment of parameter estimation in the case of CMB 
data alone uses the most general cosmology yet studied 
with both temperature and polarization information. Tak- 
ing account of differences in cosmological parameteriza- 
tions, fiducial models, and experimental specifications, we 
compare our results with past work in a series of tables. 
We generally find stronger degeneracies and hence weaker 
constraints than previous papers and propose numerical 
explanations for the discrepancies. 

We review parameter estimation methods in §0 and 
describe our parameterization of cosmology in §pl^ We 
present a way of interpreting parameter covariance in Ap- 
pendix [A| and discuss the need for careful numerical treat- 
ments in Appendix |b| In we compute the precision 
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2 



with which upcoming CMB experiments can potentially 
measure cosmological parameters. We compare our re- 
sults to previous studies in Appendix |c[ In § ^, we add 
redshift survey information and conduct several parame- 
ter studies. We explore the dependence of our results on 
our assumptions in §^|. In we discuss the necessity of 
cosmological consistency checks and highlight a number of 
possibilities. We conclude in §||. 

2. FISHER MATRIX METHODS 
Th e Fisher information matrix encodes the manner in 



which ! experimental data depends upon a set of underly- 
ing theor etical parameters that one wishes to measure (see 
Tegmark 1997b for a review). Within this set of param- 
eters, the Fisher matrix yields a lower limit to error bars 
and hence an upper limit on the information that can be 
extracted from such a data set. With the further assump- 
tions of Gaussian-distributed signal and noise, Fisher ma- 
trices can be constr ucted from the specific ations o f both 
CMB experiments ( lungman et al. 1996b; Beljak 1996b; 



Zaldarriaga i 



:;dshift 



surveys 



Selj akJ 1997 

Tegmarkl 1997a 



Kamionkowski et al.| 1997) and 
lb the extent that 



the d fita sets are independent — a very good approxima- 
tion for the cosmologies of interest — we can combine their 
constraints simply by summing their Fisher matrices. 

Suppose that the observed data are written as x\, X2, 
. . . x n , arranged as a vector x. Then suppose that the 
model parameters are p\ , P2 , ■ ■ ■ p m , arranged as a vector 
p. Let the probability of observing a set of data x given 
the true parameters p (the "fiducial model") be L(x;p). 
The Fisher matrix is then defined as 



d 2 lni 

dpidp-j 



(1) 



The Cramer-Rao inequality says that the variance of an 
unbiased estimator of a parameter pi from a data set can- 
not be less than (F~ 1 )u. In this sense, the Fisher matrix 
reveals the best possible statistical error bars achievable 
from an experiment. 

Since the error bars on our independent variables p will 
in general be correlated, they do not contain sufficient in- 
formation to calculate the errors on constructed quantities. 
A change of variables shows that the best possible error 
bar on some quantity <?(p) when marginalizing over the 
other directions that span the parameter space is 



E 



dg 
dp% 



dg 
dpj 



(2) 



More generally, the Fisher matrix transforms as a tensor 
under a change of variables in parameter space. We will 
use equation (0) to quote errors for a number of these con- 
structed quantities, such as or, the rms mass fluctuations 
in a sphere of Rh" 1 Mpc radius. 

It is at times convenient to think of the inverse of the 
Fisher matrix as a covariance matrix with an associated 
error ellipsoid. This view can be misleading. First, it rep- 
resents a degeneracy as a straight line rather than the true 
curve. For example, a CMB experiment might determine 
Vl m h 2 well but neither Q, m nor h well. The proper er- 
ror contour in the U, m -h plane would be a banana-shaped 
region along a curve of constant Q, m h 2 ; however, it will 



instead be represented as a long ellipse with the slope of 
the Sl m h 2 curve at the location of the fiducial model. 

Second, the error contours from the Fisher matrix are 
not necessarily those that would be obtained from a like- 
lihood or goodness-of-fit analysis of a particular data set. 
If the Fisher matrix errors are roughly constant across the 
error region itself (i.e. if the likelihood function is nearly 
Gaussian), then these various error bars will be compara- 
ble. This generally occurs when the error estimates are 
small compared to characteristic range over which a given 



cosmologica l parameter affects model predictions ( Zaldar- 
riaga et al. 1997a). As we will see, this is usually the 
case when CMB and redshift surveys are combined, but 
CMB data alone is subject to degenerate directions that 
surely violate the approximation. In these cases, the ques- 
tion of what confidence region to use descends into the 
murky debate between frequentists and Bayesians. If a 
strong degeneracy is present, the Fisher matrix method 
is guaranteed to find it, but different methods may dis- 
agree on the size and shape of the error region. After 
this paper was submitted, a paper investigatin g this issue 
with Mo nte Carlo methods was submitted by Efstathioi 



fc Bond (1998) 



2.1. CMB anisotropics 

Under the assumption of Gaussian perturbations and 
Gaussian noise, the Fisher matrix for CMB anisotropies 
and p olarization is (Seljak 1996b; Zaldarriaga & Seljak 
1997; |Kamionkowski et al.| 1997) 



EE 

t X,Y 



dC 



xe 



dpi 



(Cov e 



'XY 



dC 



XI' 



d Pj 



(3) 



where Cxe is the power in the I th multipole for X — T, E, 
B, and C — the temperature, iJ-channel polarization, B- 
channel polarization, and temperature-polarization cross- 
correlation, respectively. We will use Ct at times to refer to 
all the CMB power spectra together. The elements of the 
(symmetric) covariance matrix Cove between the various 
power spectra are 



(Cov^ 



TT 



-{Cre + w^Bj 2 ) 2 



(4) 



(2^+l)/ sky 

{Cove) EE = lnt {Cei + Wp^BJ 2 ) 2 , (5) 



(2^+l)/ sky ' 

(Cov,)bb - , 2 1W (C B e + wp 1 ^ 2 ) 2 , (6) 



(2^+l)/ sky ' 

(Cov^cc = — r - 7T7 [C 2 C e + {C T e + w^B^ 2 ) 



(Gove) T E 
[Cove)TC 



(2£+l)f sky 

x(C E e + w P 1 B^ 2 ) 

2 2 
(2£+l)U y Cce ' 



(7) 
(8) 



-CceiCre + w^B^ 2 ), (9) 



(2t + l)/ sky 

(Cove) EC = , 0/ ,, 2 lW Cce(C E t + w^Bf), (10) 



(Cove)TB = {Cove)EB = (Cov £ ) C s = 0. (11) 

Here B 2 is the beam window function, assumed Gaus- 
sian with B 2 = exp{-£(£ + l)6» 2 cam /81n2), where 6> boam 
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Table 1 



CMB Experimental Specifications 



Experiment 


Frequency 


^bcam 


(Tj 1 


<7p 


MAP: 


40 


28.2 


17.2 


24.4 




60 


21.0 


30.0 


42.6 




90 


12.6 


49.9 


70.7 


Planck: 


143 


8.0 


5.2 


10.8 




217 


5.5 


11.7 


24.3 



NOTES. — Frequencies in GHz. Beam size 0beam is the FWHM in 
arcminutes. Sensitivities cry and op are in fiK per FWHM beam 
(and hence Ci must be in /^K 2 ). w = (#beam CT )~ 2 is the weight 
given to that channel. 



is the full-width, half-maximum (FWHM) of the beam in 
radians. % and wp are the inverse square of the detec- 
tor noise level on a steradian patch for temperature and 
polarization, respectively. A fully-polarized detector has 
wp — 2wt- For multiple frequency channels, wB 2 is re- 
placed by the sum of this quantity for each channel. These 
formulae are derived in the / s k y = 1 case; the approxima- 
tion for / s k y < 1 including sample variance (Bcott et al 



1993) only gives the correct Fisher matrix providing that 
the power spectra have no sharp spectral features on scales 
A£ <; A# _1 , where A8 is the angular extent of the map 
in the narrowest direction (Tegmar | 1997b). Within the 
class of models we are considering, this should be an excel- 
lent approximation for MAP and Planck except at / < 3 
where sample variance is large anyway. 

We normalize th e CMB power sp ectra to COBE when 
using equation (||) ( Bunn fc White 1997). 

In Table 0, we list the experimental specifications for 
the MAP and Planck satellites used in this paper. We use 
/sky = 0.65 in all cases. We are not using the 22 GHz or 
30 GHz channels of MAP and are using only 2 of the 10 
channels of Planck. The rationale is that the statistical 
power of these channels will be used for multifrequency 
subtraction of foregrounds, leaving the full power of the 
remaining chan nels for cosmological use. We will discuss 
this further in § 3.1 



2.2. Redshift surveys 



k 2 vanish as k 
used for parameter estimation; we use 0.1/iMpc -1 as a 
default but will discuss this at length. 

Equation ( |l2] ) was derived under the approximation that 
the galaxy distribution is that of a Gaussian random field 
and that the power spectrum has no features narrower 
than the inverse scale of the survey. It also neglects edge 
effects and redshift distortions. Fortunately, the SDSS is 
both wide angle (tt steradians) and deep, so power spec- 
trum features (e.g. baryonic oscillatio ns) should be well- 



resolv ed and edge eff ects manageable (Heavens fc Taylor 



1997; [Tegmark et al.| 1998b; c.f. |Kaiser fc Peacock| 1991 
The 2dF survey, however, has a more complicated geome- 
try and will require a more careful analysis. 

We will quote our results for the Bright Red Galaxy 
(BRG) portion of the SDSS. This subsample will be in- 
trinsically red galaxies; such galaxies tend to be bright 
cluster galaxies and so the sample will reach significantly 
deeper than the primary survey. We assume the sample 
to be volume-limited with 10 5 galaxies to a depth of 1 
Gpc and to have a bias such that <7s,gai = 2. We will also 
consider the results for the main SDSS, which includes 10 6 
galaxies with a more complicated radial selection function; 
we assume this sample to have <t s ga i = 1. The value of 
08, gal affects the normalization of P(k) in equation (|l3"|); 
together with the assumed CMB normalization and spec- 
tral tilt, it implies that the fiducial model may have galaxy 
bias b ^ 1. 

We will focus entirely on the power spectrum at large 
scales, where linear theory is expected to be a good ap- 
proximation. We do this by choosing /c max to be roughly 
the s cale at w hich non-linear clustering becomes impor- 
tant (Tegmark 1997a). While data on smaller scales will 
yield very accurate measures of the power spectrum and 
higher-order correlation functions, their interpretation in 
terms of cosmological pa rame ters is m uch more compli- 
cated (Fry & Gaztanaga 1993; Peacock| 1997; Mann et al 



For the power spectrum derived from galaxy redshift 
survey s, the Fisher matrix may be approximated as QTcg-| 
mar k| 1997a) 



1998). To be conservative, we are neglecting the cosmo 
logical information on non-linear scales. SDSS should also 
reve al a wealth of in formation abou t redshift distortions 
(e.g. [Hamiltorj 1997; |Hatton fc Colej 1998) on both linear 
and non-linear scales; we will return to this in § ^. 

We assume that the galaxy bias is linear on these large 
scales. This has so me t heoretical j ustific ation (|Colcs| 1993 ; 
|Fry fc Gaztanaga 1993; Weinberg 1995; Bcherrer & Wein- 



dlnP(fc) dlnPfk) , ,k 2 dk 
V c s{k)- 



V eS (k) 



dpi 
n(r)P{k) 
1 + n(r)P(k) 



d Pj 



(2nf 



(12) 
(13) 



where n(r) is the survey selection function, i.e. the ex- 
pected number density of galaxies at the location r. P(k) 
is the model galaxy power spectrum. V c s{k) is the effec- 
tive volume of the survey, properly weighing the effects 
of shot noise in undersampled regions. fc m i n is the mini- 
mum wavenumber to which the survey is sensitive, but in 
practice the numerical results are virtually unchanged by 
taking fc m ; n = since both V c s and the phase space factor 



berg 1997; Mann et al. 1998). More important, this as- 
sumption will be stringently tested by the SDSS and other 
large surveys. Galaxies of different morphologies or type 
will have different levels of bias, but linear bias predicts 
that the ratios of the various power spectra should be con- 
stant on large scales. Scale-dep endent bias has been de- 
tected on small scales (Peacock 1997), but this does not 
test the linear bias assumption we are making here. In 
addition, redshift distortions may be able to probe any 
scale-dependence of bias on large scales. 

3. PARAMETERIZED COSMOLOGY 

We adopt a 13-dimensional parameterization of the adi- 
abatic CDM model. Our independent variables include 
the matter density £l m h 2 , the baryon density f2s h 2 , the 
massive neutrino density f2„/i 2 , the cosmological constant 
Q\, and a curvature contribution £Ik- Here, the Hubble 
constant is written as H = lOOh km s _1 Mpc -1 . These 



4 




/ / 

Fig. 1. — Our most used model along with 1— a band-power error bars from MAP and Planck. The model is Q m = 0.35, h = 0.65, Qg = 0.05, 
Q A = 0.65, ft„ = 0.0175, t = 0.05, n s (fc fld ) = 1, and T/S = a = 0. a) The temperature power spectrum AT = [£(£ + l)C T e/^] 1/2 - b) 
The E-channel polarization power spectrum AT = [t(t + 1)Ce£ /27T] 1 / 2 ; the panel shows a blow-up of the large-angle feature caused by 
reionization. MAP errors are the lighter, larger boxes; Planck errors are the darker, smaller boxes. The bands reflect an averaging over many 
l\ the actual experiments will have finer I resolution (and correspondingly larger errors). MAP will be able to average the polarization bands 
together to get a marginal detection at I fa 150, but Planck can trace out the full curve. Note that because polarization and temperature are 
correlated, the significance of detecting a change in parameters using both data sets is not simply given by the combination of errors from 
each. 



definitions imply that the total matter density in units of 
the critical density £l m = 1 — fl\ — £Ik, that the Hub- 
ble constant h = yj (n m h 2 )/tt m , and that the CDM den- 
sity Ocdm^ 2 = Q m h 2 — ilsh 2 — VL v h 2 . We use a sin- 
gle species of massive neutrinos, so the neutrino mass is 
m v « 94fl„h 2 eV (we take units in which the speed of light 
is unity). 

We include an unknown optical depth r to reionization, 
implemented as a rapid and complete ionization event at 
the appropriate (small) redshift. We also allow the pri- 
mordial helium fraction to vary but assume that the infor- 
mation from direct abundance measurements can b e rcp- 
resented by a Gaussian prior Y v = 0.24±0.02 (1-er) ( Bond| 
et al.[]l997; |Schramm fc Turncr| 1998). 

We use an initial power spectrum of the form 



= A 2 (fc/fcfid)" s(fe) - 4 (14) 



with 



n s (k) = n s (/cfid) + arn(fc/fc fid ) (15) 

for the fluctuations in the gravitational potential. Here, a 
is a logarithmic running of the tilt around a fiducial scale 
fcfid = 0.025 Mpc -1 . The density power spectrum is equal 
to the potential power spectrum times (fc 2 + SQkHq ) 2 ; 
the simplest open inflationary models predict a power- 
law in the potential, not in the density. In a flat uni- 
verse, the initial density power spectrum takes on the usual 
A 2 s {k/kf i< \) ns form. The present-day density power spec- 
trum P(fc) of course differs by the square of the transfer 
function. 

Note that equivalent parameters for a different choice 
of fcfid could be mapped into the parameters of equation 
(jl4]); hence the value of fc nd is immaterial. However, be- 
cause of the running of the tilt, the value and hence the 
error bars on the tilt itself become scale-dependent. We 



will quote values at both the Hubble wavenumber and at 
fcfid- The errors on a are scale- independent for the param- 
eterization of equation ([u]) . As shown in Appendix [XJ for 
a given experiment and fiducial model, there is a "pivot" 
wavenumber fc P i VO t for which the errors on ns(fc p i vo t) and 
a are uncorrelated and the errors on ns(^pivot) are equal 
to the errors on ns when a is held fixed. As one might 
expect, fcpivot fohs near the center of the observable range 
of wavenumbers, generally not too far from our choice of 
fcfid- At other wavenumbers, the uncertainties on ns(fc) 
are larger and correlated with a. 

We allow tensor perturbations with a normalization T/S 
equal to the ratio of the I — 2 temperature anisotropics of 
the tensors and scalars. Note that this is proportional but 
not equal to the ratio A\j A 2 S that enters into the inflation- 
ary constraints. In particular, cosmological parameters en- 
ter i nto T/ S due to the evolution of quadrupole anisotro- 
pics (Knox 1995; Turner fc White 1996). We parameterize 



the power-law exponent of the tensor input spectrum as 
tit] this allows us to probe whether these data sets can 
test the inflationary consistency relation. However, for 
fiducial models with T/S = 0, excursions in tit occur at 
T/S = 0, thereby yielding a zero derivative. Hence, in this 
limit, tit is not a physically meaningful parameter and our 
parameter space is effectively reduced to 12 dimensions. 

Finally, we allow the scalar normalization to vary and 
include an unkno wn line ar bias. We describe our normal- 
ization choice in § B.2.2 ; this choice can affect individual 
derivatives but does not affect marginalized errors. The 
linear bias b is defined by P ga \ = b 2 P mass . We also quote 
results for (3 = f2^ 6 /b to fac ilitate c omparisons to results 



from peculiar velocity data ( Peebles! 1980). 

We discuss numerical issues involved with constructing 
derivatives with respect to these parameters in Appendix 
|b[ Here we simply note that our results in Table || are 
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lar diameter distance degeneracy and the reionization 
normalization degeneracy. Although these were discussed 
in previo us papers, none o f the ta bles in Jungman ct al 



k (h Mpc- 



Fig. 2. — The power spectrum for the model of Figure hi. SDSS 
BRG 1— <r error bars are superposed. Dashed boxes are for k > 
fcmax = O.lh Mpc -1 ; we neglect information from these scales unless 
otherwise noted. 



stable to 10% when step-sizes are halved and that we have 
spot-checked elsewhere with similar results. Stability im- 
proves as parameter degeneracies are removed, either by 
the addition of other data sets or by external priors. 

3.1. Fiducial Models 

The Fisher matrix formalism asks how well an exper- 
iment can distinguish the true ("fiducial") model of the 
universe from other models. The results clearly depend 
upon the fiducial model itself. We will quote results for 
a number of models, but most of our studies are based 
around a low-density, geometrically flat ACDM model 
with n m = 0.35, h = 0.65, n B = 0.05, n A = 0.65, 
n„ = 0.0175 K = 0.7 eV), r = 0.05, n s (fc fid ) = 1, 
and T/S = 0, in the region of parameter space favored by 
recent data. We display the CMB and matter power spec- 
tra for this model in Figures [j] and |^. Note that although 
our fiducial model is flat, we do not restrict our excursions 
to flat models: we vary all 12 parameters (including £Ik) 
simultaneously unless explicitly stating otherwise. Numer- 
ical issues are the reason for the small but non-zero values 
of 0„ and t here (see Appendix 0). We will alter the 
fiducial values o f f2 #, f2„, and T/S in various parame- 
ter studies. 



In 



5.3 



we will consider an Q, m = 1 SCDM 
model, a Q m = 0.35 open CDM model, and a f2 m = 0.2 
flat model, a = and nx = in all models. 

4. CMB ALONE 

We begin by presenting the constraints that could in 
principle be achieved by CMB satellites without galaxy 
power spectrum information. Table shows the results 
on a variety of quantities using a fl m = 0.35, fl A = 0.65 
CDM model as our fiducial model. As expected, CMB 
data alone provides excellent constraints on il^/i 2 , fl m h 2 , 
ns(fcfid) and a, as well as on a combination of fix and 17a- 
With polarization information, it also strongly constrains 
the reionization optical depth r. 

However, the remaining parameters are nowhere near 
"percent level" in accuracy, even for Planck. This is 
because of various degeneracies, particularly the angu- 



Copeland et al 



(1996b), Zaldarriaga ct al. (1997), |Bond et al. (1997) or 



1998) include the case where SIk and Q, A 
are measured simultaneously. As is well known, a com- 
bination of changes in Oa and VLk at fixed $l m h 2 and 
fish 2 can keep the angular location and morphology of 
the acoustic peaks fixed. This ambiguity ca uses th e value 
of h, f2 m , and £l m h to be poorly constrained (EHT); more- 
over, these uncertainties propagate into the rms density 
fluctuations <jr due to the redshift-space definition of the 
tophat radius R (h^ 1 Mpc). Any additional parameter 
that affects the angular distance relation, e.g. variations 
in the equation of state of the missing energy, creates a 
similar degeneracy. 

Adding polarization information to the CMB helps for 
most quantities. Much of the improvement comes from 
the ability of polarization infor mation to isolate re ioniza- 
tion and tensor contributions (Zaldarriaga et al. 1997), 



thereby separating the large-angle temperature effects of 
Q A and VLk- Table ^ illustrates that removing this ambi- 
guity substantially reduce errors on other quantities even 
with the polarization sensitivity of MAP. 

Because the temperature-polarization correlation is not 
complete, polarization also provides independent informa- 
tion on the high-redshift, degree-scale acoustic oscillations 
to combat sample variance. This is best illustrated by fix- 
ing r, T/S, and Qk, thereby eliminating the need to em- 
ploy the large-angle polarization signal to break degenera- 
cies. For the model in Table |^, adding polarization infor- 
mation to the temperature data in this restricted parame- 
ter set improves errors on Q m h 2 , fish 2 , fl A , and ns(fcfid) 
by about 40%-70% for Planck. The improvement is only 
8% for MAP, so polarization maps of this sensitivity are 
useful only through their constraints on r and T/S as dis- 
cussed above. 

Parameter estimation with the CMB alone has been dis- 
cussed in the past in smaller parameter spaces and/or 
without polarization information. As discussed in Ap- 
pendix |C, our results disagree with several studies in the 
literature, generally by giving larger error bars. 

5. COMPLEMENTARITY 
5.1. Adding Redshift Survey Data 

The addition of information on the matter power spec- 
trum at the precision available within the SDSS can make 
a substantial improvement in the error bars on key cos- 
mological quantities. As seen in Table |, one gets large 
improvements on h, O m , and the relatedquantities f2 m /i, 
Qaj and and moderate improvements in other quan- 
tities, particularly VL v h 2 . Of course, as one improves the 
quality of the CMB data set, the fixed level of SDSS input 
gets less and less important. 

The most striking improvement allowed by measure- 
ment of the matter power spectrum is the breaking of 
the angular diameter distance degeneracy. Since Q\ and 
VLk shift the acoustic peaks in opposite directions, to 
trade one off the other requires substantial changes in 
Q m = 1 — tt\ — fix. Because this variation must be done 
at fixed fl m h 2 to maintain the peak morphology, h varies 
strongly in the degenerate direction. A measurement of h 
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Table 2 



Marginalized Errors for ACDM for various CMB experiments. 







CMB alone 




SDSS 




CMB- 


hSDSS 






MAP 


Planck 


clloilC 


MAP 


Planck 


Quantity 


Temp 


1 +f 


Temp 


1 +f 




Temp 




Temp 


■ | ■ 1 r> 
1 +f 


h 


1.3 


0.22 


1.1 


0.13 


1.3 


0.030 


0.029 


0.025 


0.022 




1.4 


0.24 


1.2 


0.14 


0.23 


0.042 


0.036 


0.035 


0.027 




0.47 


0.078 


0.40 


0.046 


0.59 


0.024 


0.018 


0.015 


0.010 


n A 


1.1 


0.19 


0.96 


0.11 


CO 


0.056 


0.042 


0.036 


0.024 




0.31 


0.055 


0.26 


0.030 


CO 


0.022 


0.015 


0.007 


0.005 


ln(C m /i 2 ) 


0.20 


0.095 


0.064 


0.018 


4.5 


0.11 


0.077 


0.040 


0.016 


ln(n B h 2 ) 


0.13 


0.060 


0.035 


0.010 


5.6 


0.074 


0.050 


0.026 


0.010 


m v (eV) oc S!j,/i 2 


0.89 


0.58 


0.58 


0.26 


9.1 


0.38 


0.33 


0.31 


0.21 


*p 


0.020 


0.020 


0.018 


0.013 




0.020 


0.020 


0.017 


0.013 




0.11 


0.048 


0.041 


0.008 


1.1 


0.064 


0.040 


0.028 


0.008 




0.32 


0.17 


0.18 


0.039 


4.1 


0.23 


0.14 


0.13 


0.038 


a 


0.030 


0.018 


0.015 


0.004 


0.42 


0.022 


0.015 


0.011 


0.004 


lnP s (fc fld ) = In A| 


1.4 


0.43 


1.1 


0.073 


CO 


0.61 


0.36 


0.36 


0.069 


In P$(i? ) 


1.8 


0.71 


1.3 


0.16 


CO 


0.91 


0.61 


0.55 


0.15 


T/S 


0.48 


0.18 


0.35 


0.012 




0.24 


0.16 


0.15 


0.012 


r 


0.69 


0.022 


0.59 


0.004 




0.27 


0.021 


0.21 


0.004 


In o-g 


0.48 


0.14 


0.42 


0.057 


CO 


0.27 


0.070 


0.22 


0.044 


ln(o"5o/o-8) 


0.86 


0.15 


0.75 


0.093 


0.27 


0.029 


0.028 


0.024 


0.020 


ln/3 










CO 


0.27 


0.068 


0.20 


0.027 


In 6 










CO 


0.28 


0.087 


0.23 


0.062 



NOTES.— n m = 0.35, tt B = 0.05, H A = 0.65, h = 0.65, T/S = 0. n T = and cannot vary. All errors arc I - a. k maK = O.lhMpc -1 . 
Blank entries indicate that the parameter doesn't affect the observables for the data set. Infinite entries indicate that the parameter affects 
observables but is not constrained due to degeneracies; in particular, the growth factor, normalization, and bias are all degenerate. 



thereby breaks the degeneracy. 

The combination of CMB and galaxy survey data can 
do just that. Once fl m h 2 and Ugh 2 are well-determined 
by the CMB acoustic peak morphology, the real-space 
power spectrum is known. Varying h causes this pattern 
to move in redshift space. Hence, the more features that 
are present in the matter power spectrum, the more accu- 
rately one can measure h and so break the degeneracy. 

Within the types of CDM models considered here, the 
best source of power spectrum features are the o scillat ions 
i mpr essed by baryonic oscillations. As shown in EHT and 
§ even a 10% baryon fraction causes large enough fea- 
tures for SDSS to clamp down on h and Q m . Indeed, the 
resulting error bars are nearly as good as if the universe 
were assumed flat in many cases, thereby breaking the de- 
generacy by fiat. 

In other sectors, the gains are more modest. Improve- 
ments approaching a factor of 2 are possible in fl m h 2 and 
ns (km), particularly for MAP without polarization. Po- 
larization tends to allow the CMB to dominate the con- 
straints, but this does depend on the value of fc ma x and on 
being able to extract low-£ cosmological information. An 
important complementary aspect of the data sets enables 
the de termination of the mass of cosm olog ical neutrinos 



(|HETJ). We will discuss this further in §p.4. 

As shown in Table g, SDSS power spectrum information 
alone flounders in this large parameter space. Details of 
the shape do depend on cosmology, but the features are 
not well enough detected for this low baryon fraction and 
small fc max (see Figure ||) . 

In all cases, the improvements depe nd s omewhat on the 
value of fc max , which we will vary in § 3.3, Reverting from 



the deeper BRG survey to the SDSS main survey degrades 
the performance on h and related quantities by about 30%. 

5.2. The Role of Baryons 

Table || shows the results as a function of the baryon 
fraction. As the baryon fraction increases, the addition 
of SDSS information becomes more and more helpful, h 
and the related quantities f2 m , VIk, and as are the 
most affected, but even traditional CMB quantities such as 
Cl m h or flgh 2 see marked improvement at high flg/Q m . 

The driving physical effect behind these gains is the 
structure that develops in the matter power spectrum due 
to the high-redshift acoustic o scillations imp rinted by a 
non-negligible bary on fr action (Peebles fc Yvj 1970; Sun-| 
lyaev fc ZePdovichl 1970 ; |Holtzman| 1989; [Hu fc Sugiyama) 
1996; jEiscnstcin L Hu| 1998a). In CDM cosmologies, the 



baryons and photons oscillate on sub-horizon scales prior 
to recombination, while the CDM perturbations simply 
grow. The event of recombination catches the oscillations 
at various phases and creates the acoustic peaks we see 
in the CMB spectrum. The perturbations in the bary- 
ons also share this oscillatory history, but in trace-baryon 
cosmologies it is erased as the baryons fall into the more- 
evolved CDM perturbations. When the baryon fraction 
is non-negligible, however, the equilibration of the baryon 
and CDM perturbations is not completely one-sided, and 
the final power spectrum P(k) retains an imprint of the 
acoustic oscillations. The resulting morphology consists of 
a sharp break in the power spectrum followed by a damped 
series of wiggles; the whole pattern has a characteristic 
scale, known as the sound horizon, which is the distance a 
sound wave could travel prior to recombination. 
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Table 3 



Marginalized Errors as Function of fl B . 





fig 


= 0.005 


n B 


= 0.02 


n B 


= 0.05 


n B = 


0.10 


Quantity 


MAP 


+SDSS 


MAP 


+SDSS 


MAP 


+SDSS 


MAP 


+SDSS 


h 


0.34 


0.12 


0.27 


0.091 


0.22 


0.029 


0.23 


0.013 


Q-rri 


0.36 


0.12 


0.29 


0.086 


0.24 


0.036 


0.25 


0.018 


k 

1 Lm >v 


0.12 


0.039 


0.093 


0.029 


0.078 


0.018 


0.084 


0.009 




0.28 


0.098 


0.23 


0.070 


0.19 


0.042 


0.20 


0.022 


n K 


0.082 


0.029 


0.065 


0.027 


0.055 


0.015 


0.056 


0.008 


\n(n m h 2 ) 


0.091 


0.089 


0.11 


0.101 


0.095 


0.077 


0.073 


0.034 


lll\\L B ll J 


U.UOo 


U.UD I 


U-UOl 




U.UDU 


u.uou 


U.UDZ 




m v (eV) oc Qvh 2 


0.60 


0.55 


0.60 


0.46 


0.58 


0.33 


0.74 


0.23 


ns(*fld) 


0.046 


0.042 


0.031 


0.030 


0.048 


0.040 


0.055 


0.027 


a 


0.032 


0.020 


0.023 


0.017 


0.018 


0.015 


0.026 


0.015 


lnP 4 ,(fc fld ) = In A| 


0.39 


0.35 


0.30 


0.28 


0.43 


0.36 


0.47 


0.23 


T/S 1 


0.11 


0.11 


0.15 


0.14 


0.18 


0.16 


0.17 


0.13 


T 


0.036 


0.033 


0.026 


0.024 


0.022 


0.021 


0.022 


0.020 


In as 


0.19 


0.17 


0.20 


0.14 


0.14 


0.070 


0.14 


0.046 


ln(o-5o/eT 8 ) 


0.27 


0.031 


0.19 


0.026 


0.15 


0.028 


0.13 


0.026 


In /3 




0.074 




0.069 




0.068 




0.050 


In 6 




0.31 




0.23 




0.087 




0.046 



NOTES. — All models have f2 m = 0.35, Oa = 0.65, h = 0.65, and T/S = 0. n T = and cannot vary. All errors arc 1 - a. CMB data is for 
MAP with temperature and polarization information. SDSS column uses information up to fc max = 0.1/iMpc -1 as well as CMB data. 



As discussed in EHT , once CMB data yields £l m h 2 and 
ilsh 2 , the physical size of the sound horizon is known. 
Measuring this scale in a redshift survey thereby allows 
a comparison that gives the Hubble constant. Once the 
baryon fraction is about 10% or greater, SDSS can detect 
the baryonic features and thereby collapse the error bars 
on h and related quantities. 

In Table |3|, the case of negligible baryons = 0.005 al- 
lows SDSS to make modest improvements over MAP alone. 
The factor of 3 improvement on h and fl m comes from the 
two scales that remain as £Ib gets small: the scale of the 
horizon at matter-radiation equality, which is proportional 
to £l m h, and the scale of the horizon when the massive 
neutrinos become non-relativistic. The latter is a result of 
using a 5% neutrino fraction; the former is always present 
but is inaccurately measured due to confusion with spec- 
tral tilt. 

Baryon fractions exceeding 10% a re strongly f avored 
by observatio ns of cluster X-ray gas ( White et alj 1993a; 
David et al.] 1995; |White fc Fabiarj 1995; |Evrard| 1997). 



Lyman a forest theories also favor high baryon densities 
( Weinberg et al. 1997); when combined with the ge neral 



obser vatio nal pr eference for low matter densities (e.g. Bah . 
call e1_aL 1997; |Carlberg et al. 1997ab), this yields a high 
baryon fraction. Hence, it seems likely that baryonic fea- 
tures will be prominent enough in the matter power spec- 
trum t hat SPSS and perhaps 2dF will be able to detect 
them ( [Tegmark| 1997a; |Goldberg fc Strauss| 1998). 

If the baryon fraction is yet higher, perhaps 20%, then 
the detailed morphology of the baryon features can be 
studied well enough to tighten constraints on other, non 
/i-related, parameters. The final columns in Table || show 
that a baryon fraction of 28% could allow SDSS to make 
a factor of 2 improvement in fl m h 2 , Sigh 2 , and JT-s(fcfid) 
over MAP with polarization. 



5.3. Variations in Q m and D,k 

We show the marginalized errors for three other choices 
of VL m and VLk in Table [l| The results are similar to those 
of the previous sections, with some expected variations. 
Performance with SDSS on h and fl m is worse (better) 
in the Q m = 1 SCDM (f2 ro = 0.2 ACDM) model because 
the baryonic oscillations in the matter power spectrum 
get horizontally shifted to less (more) favorable locations 
relative to the fixed cutoff scale of fc max = O.lh Mpc -1 . On 
the other hand, the degree of non-linearity at that scale is 
not fixed. Abundances of rich clusters require that low-f2 m 
cosmologies h ave a higher erg than high- f2 m cosmo logies 
(IWhite et al.| 1993b; |Viana fc Lid'dTq' 1996; [Eke et al.| 1997; 
Pen 1997); this suggests that a given level of non-linearity 
should occur at high values of k in higher Q, m cases. In 
other words, had we fixed /c max by requiring that cluster- 
normalized fluctuations reach a particular amplitude on 
that scale, the number of baryon oscillations in the linear 
regime, and hence the performance of SDSS, would have 
remained more constant. 

As expected, the fl m = 0.35 open cosmology does show 
somewhat worse performance than its flat cousin for MAP 
data, presumably because the acoustic peaks have been 
shifted to smaller scales, leaving less structure to be re- 
solved by the beam. This effect is smaller for Planck — 
while fewer peaks are detected, the sample- variance errors 
on scales around the sound horizon are improved. 

5.4. Massive Neutrinos 

Massive neutrinos present particular problems for Fisher 
matrix analyses of these data sets. Varying the neutrino 
mass changes the free-streaming scale, below which their 
rms velocity prevents them from clustering. As Q u (or 
equivalently m u ) approaches zero, the neutrinos become 
fully relativistic and the free-strea ming scale of the neu- 
trinos approaches the horizon scale (Bond fc Szalay 1983; 
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Table 4 



Marginalized Errors for Various Cosmological Models. 







SCDM (0 


-ni) 

B — "* / 






OCDM (Q 


m = 0.35) 




ACDM (Q m = 0.2) 






CMB alone 


CMB4 


SDSS 


CMB alone 


CMB- 


r-SDSS 


CMB alone 


CMB+SDSS 


Quantity 


M A P 


Pla,nck 


MAP 


Planck 


A/TAP 


Planck 


1\/T A P 


Planck 


M A P 


Planck 


A/TAP 


Planck 


h 


0.39 


0.16 


0.054 


0.036 


0.24 


0.11 


0.038 


0.020 


0.30 


0.14 


0.028 


0.023 




1.5 


0.64 


0.20 


0.15 


0.24 


0.12 


0.038 


0.023 


0.15 


0.074 


0.015 


0.012 


O h 


0.38 


0.16 


0.055 


0.038 


0.075 


0.038 


0.020 


0.008 


0.062 


0.030 


0.011 


0.005 




1.3 


0.55 


0.18 


0.13 


0.49 


0.24 


0.11 


0.051 


0.12 


0.055 


0.021 


0.010 


Q K 


0.22 


0.090 


0.046 


0.020 


0.25 


0.13 


0.072 


0.028 


0.039 


0.019 


0.010 


0.003 


ln(C m /i 2 ) 


0.082 


0.014 


0.076 


0.013 


0.14 


0.010 


0.103 


0.010 


0.097 


0.020 


0.078 


0.015 


ln(Q s /i 2 ) 


0.044 


0.009 


0.043 


0.008 


0.101 


0.008 


0.053 


0.007 


0.062 


0.011 


0.050 


0.011 


mi/ (eV) tx Q v h 2 


1.4 


0.35 


0.68 


0.31 


0.82 


0.15 


0.38 


0.14 


0.59 


0.25 


0.22 


0.15 


«s(%d) 


0.042 


0.008 


0.037 


0.007 


0.075 


0.006 


0.038 


0.006 


0.047 


0.009 


0.038 


0.008 


o 


0.016 


0.004 


0.015 


0.004 


0.043 


0.003 


0.022 


0.003 


0.027 


0.004 


0.016 


0.004 


lnP s (fc fld ) = In A 2 


0.38 


0.067 


0.34 


0.066 


0.59 


0.050 


0.33 


0.049 


0.42 


0.077 


0.34 


0.067 


T/S 


0.24 


0.019 


0.22 


0.019 


0.29 


0.016 


0.25 


0.016 


0.13 


0.009 


0.12 


0.009 


T 


0.020 


0.004 


0.020 


0.004 


0.042 


0.005 


0.041 


0.005 


0.025 


0.005 


0.022 


0.004 


Intra 


0.34 


0.14 


0.101 


0.050 


0.31 


0.052 


0.090 


0.034 


0.17 


0.061 


0.052 


0.036 


ln(o-5o/o-g) 


0.42 


0.17 


0.042 


0.034 


0.23 


0.071 


0.037 


0.016 


0.13 


0.076 


0.027 


0.021 


ln/3 






0.051 


0.031 






0.074 


0.025 






0.076 


0.034 


In 6 






0.18 


0.12 






0.11 


0.052 






0.063 


0.047 



NOTES.— SCDM model has fl m = 1, fi B = 0.10, H A = 0, h = 0.5, and T/S = 0. OCDM model has Q, m = 0.35, f2 s = 0.05, H A = 0, h = 0.65, 
n s = 1.25, and T/S = 0. ACDM model has fi m = 0.2, tt B = 0.03, C A = 0.8, h = 0.8, and T/S = 0. n T = and cannot vary. All errors are 
1 — it. fc max = 0.1/iMpc -1 . CMB experiments include temperature and polarization information. 



pairs of models that bracket the mass transition for some 
scale of interest? Second, since the derivative is changing 
significantly even for small values of £l v h 2 , how should we 
interpret the resulting error bars on Q, v h 2 l The difficulty 
is that the log- likelihood function is not well-approximated 
by a quadratic expansion around — relative to the 
quality of the constraints. 

This situation is shown in Figure |^. The temperature 
C'e derivative with respect to the massive neutrino density 
is a strong function of the fiducial value of Sl v . For small 
Sl„, the derivative has a large peak at £ s=s 150, but for 
larger the feature goes away entirely. This means that 
if one began at f2„ = and increased the neutrino mass, 
the power at £ = 150 would change rapidly at first and 
then saturate. A Fisher matrix analysis around Cl v ~ 
would show strong limits on Q„, while an analysis around 
Sl v /VL m « 0.2 would show weak limits. 

The derivative of the matter power spectrum with re- 
spect to the neutrino mass shows a similar behavior: it is 
non-zero only below a break scale that shifts to small k 
as Vl v h 2 gets smaller. Because of degeneracies, the break 
itself needs to be detected. Redshift-survey data is less 
restrictive at small k. Hence, whereas CMB data is more 
restrictive at low Q v , redshift survey data is more restric- 
tive at high f2„. 

We are interested in what neutrino masses can be dis- 
tinguished from zero. Because of the above trend, even if 
a data set preferred Sl u = 0, the upper limit on the neu- 
trino mass would be underestimated. To be conservative, 
we run at =/= instead; the significance at which f2„ = 
can be excluded by the data set is then underestimated. 

In Table ^| we present marginalized errors as a function 
of Q v within a f2 m = 1 CDM model. MAP data alone does 
not detect the mass of the neutrino in any of the cases. 
In particular, as Q, u increases, the angular scale at which 




50 100 



500 



1 



Fig. 3. — The derivative d(ln Cti)/ d£l u h 2 as a function of neutrino 
fraction Q v /Q m - Top to bottom (for the first peak): 0.5%, 1%, 2%, 
5%, and 20%. The cosmology is the U m = 0.35 ACDM model, but 
with t = 0.1, for a variety of f!„. 



Holtzman 1989; Ma fc Bertschinger 1995; Dodclson ct al_ 
1996; |Hu fc Eiscnstcirj 1998). This means that at large 



scales, it takes only a small upward variation of the neu- 
trino mass from zero to bring the neutrinos out of the 
free-streaming regime. As the power spectra phenomenol- 
ogy differs between these two physical regimes, we expect 
that the derivative of P(k) or C; with respect to il„h 2 
will be sharply different above and below this small mass 
threshold, which is itself a function of scale. 

This raises both a practical and a theoretical concern. 
First, given numerical noise, is it possible to calculate the 
derivative as fl u goes to zero, or must one always difference 
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Table 5 



Marginalized Errors as a Function of fi„ 







Q 


„ = 0.01 






n„ = 0.05 




n„ = 0.10 


n„ 


= 0.20 








\7Tlv 


= 0.24 e 


V) 


( 


<7l„ = 1.2 


eV) 


r 


ti„ = 2.4 


eV) 


{vriu '- 


= 4.7 eV) 






CMB 


+SDSS 


CMB 


+SDSS 


CMB 


+SDSS 


CMB 


+SDSS 


Quantity 




T+P 


0.1 


0.2 


T+P 


0.1 


0.2 


T+P 


0.1 


0.2 


T+P 


0.1 


0.2 


m v (eV) oc 


VL v h 2 


1.1 


0.81 


0.42 


1.4 


0.68 


0.30 


5.5 


0.84 


0.31 


6.8 


1.7 


0.38 


ln(Q m /i 2 ) 




0.11 


0.101 


0.082 


0.082 


0.076 


0.063 


0.067 


0.065 


0.053 


0.084 


0.059 


0.044 


h 




0.38 


0.065 


0.020 


0.39 


0.054 


0.018 


0.39 


0.046 


0.017 


0.38 


0.039 


0.015 






1.5 


0.28 


0.092 


1.5 


0.20 


0.068 


1.5 


0.16 


0.055 


1.5 


0.12 


0.039 


m,, (eV) oc 


VL v h? 


0.53 


0.47 


0.29 


0.35 


0.31 


0.19 


0.30 


0.28 


0.17 


0.31 


0.31 


0.17 


\n(n m h 2 ) 




0.028 


0.025 


0.018 


0.014 


0.013 


0.011 


0.010 


0.010 


0.010 


0.009 


0.009 


0.009 


h 




0.16 


0.048 


0.016 


0.16 


0.036 


0.014 


0.16 


0.028 


0.012 


0.16 


0.021 


0.009 






0.64 


0.21 


0.075 


0.64 


0.15 


0.058 


0.64 


0.11 


0.048 


0.64 


0.081 


0.034 



NOTES. — Qm = 1, f^s = 0.1, Q\ = 0, h = 0.5, and T/S = for all columns. We assume one species of massive neutrino, tit = and does 
not vary. All errors are 1 — a. All columns include CMB information on temperature and polarization; the top four lines are for MAP, the 
bottom four for Planck. SDSS columns include information to fc max = 0.1/iMpc or /c max = 0.2ft Mpc" 1 , as noted. 



Table 6 



Marginalized Errors as a function of Tensor-to-Scalar Ratio 









T/S = 1.0 






T/S = 0.3 


T/S = 0.1 




MAP 




Planck 




Planck 


Planck 


Quantity 


T+P 


SDSS 


no B 


T+P 


SDSS 


SDSS 


SDSS 


T/S 


2.0 


0.79 


0.45 


0.42 


0.28 


0.13 


0.064 


X = T/S + xn T 


0.65 


0.54 


0.20 


0.16 


0.14 


0.055 


0.032 


Tlx 


0.40 


0.16 


0.12 


0.100 


0.080 


0.14 


0.24 


™sOfid) 


0.065 


0.047 


0.012 


0.012 


0.011 


0.009 


0.008 


ns(Ho) 


0.45 


0.21 


0.059 


0.055 


0.051 


0.043 


0.040 


a 


0.048 


0.022 


0.005 


0.005 


0.005 


0.004 


0.004 


lnP s (fc fld ) = lnA| 


0.59 


0.43 


0.11 


0.103 


0.095 


0.080 


0.074 


lnP$(tf ) 


1.4 


0.81 


0.27 


0.25 


0.23 


0.18 


0.17 


Value of x 


4.9 


3.6 


3.4 


3.7 


3.1 


0.85 


0.23 



NOTES. — Um = 0.35, Q B = 0.05, Q A = 0.65, h = 0.65, n T = 0. All errors are 1 - a. fc max = Ori/iMpc" 1 . SDSS columns include CMB data 
with polarization including B channel. The "No B" column presents results in which B-channel polarization has been ignored. The B channel 
is of negligible importance for MAP. x has been chosen so as to minimize the error on the quantity X; for this choice, X is uncorrelated with 
riT and the error is equal to the error on T/S if n-p were held fixed. Without CMB polarization information, errors on T/S from the MAP or 
Planck are much larger, but errors when combined with SDSS are similar to those in the MAP+SDSS column. 



the neutrinos have an observable signature drops below 
the resolution of the experiment. Planck can continue to 
track the effect and therefore can detect the neutrino mass. 
In both data sets, as fi„ increases, £l m h 2 errors tend to 
decrease and errors on the tilt sector remain unchanged. 

With SDSS data added to MAP, n eutrin o masses ex- 
ceeding ~ 1 eV would be detectable ( |HET| ). SDSS also 
improves the limits of Planck by up to a factor of 2. Sen- 
sitivity to lower masses is possible for O m < 1. Unlike 
the situation with the ACDM models presented elsewhere 
in this paper, fc max ~ 0.2/iMpc _1 may be appropriate for 
the neutrino signatures of SCDM: these models have lower 
normalization (e.g. cluster abundances suggest erg ss 0.5 
not 1.0) and thus less non-linear contamination. Also, be- 
cause the neutrino signature is not oscillatory, it may be 
harder for nonlinear dynamics to wash it out. 

As £l v increases, performance on Q m and h from the 
combination of SDSS and CMB also improves. Like the 
baryons, massive neutrinos impress a scale on the mat- 



ter power spectrum: there is a break in the spectral index 
around the scale of the horizon at the epoch when the neu- 
trinos become non-relativistic. This physical scale depends 
upon fl m h 2 and Q„/i7 TO ; once these two are determined, 
one gains leverage on h and f2 m by detecting the scale in 
rcdshift space. 

We have run f2„/J7 TO = 0.05 unless otherwise specified. 
With CMB and redshift survey data, this model is distin- 
guished from Q„ = at 2 — 4a. Larger neutrino fractions 
rely more and more on the matter power spectrum data 
to provide constraints. 

5.5. Tensors 

Interpreting the errors on the tensor-to-scalar ratio T/S 
requires special care since we have at present no guide 
as to its value in the real univ erse. Althou gh in power- 
law inflation T/S ~ 7(1 — ns) ( Davis ct al. 1992), many 
inflationary models predict T/S ~ (Lyth 1997) and in 
general slow-roll inflationary models obey T/S « —IriT 
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(see e.g. |Liddle fc Lyt"h| 1993). For this reason, we conduct 
a parameter study of T/S here. 

T/S Sequence — As T/S increases in the fiducial model, 
the errors on most quantities degrade slightly due to the 
reduction of small-angle anisotropy signal at fixed COBE- 
normalization. With information from SDSS, the degra- 
dation is smaller yet. However, the errors on T/S and 
tit are strong functions of T/S. Errors on nj- and frac- 
tional errors on T/S decrease as T/S increases, the latter 
of course becoming infinite as T/S — > 0. Including SDSS 
information significantly reduces the error bars even in the 
case of Planck; for example, better controlling the value 
of J1a and fix substantially helps Planck to sort out the 
various larg e-angl e signals and ap proach t he fixed cosmol- 
ogy limit of Knox (1995) (see also Kinney| 1998). We show 
results for T/S = 1, 0.3, and 0.1 in Table |6[ Remember, 
however, that each of these assumes nr = in the fidu- 
cial model; the inflationary consistency relation predicts 
nx < for T/S 0, in which case the signal would be 
slightly overestimated. 

T/S^ut plane — In the presence of a varying ny, the 
value and hence the error bars on the tensor-to-scalar ratio 
become scale-dependent. The ratio of quadrupole power 
T/S is a pessimistic choice because the ability of the data 
sets to constrain tensor contrib utions comes at somewhat 
smaller angles (£ « 50, see also White 1996). Hence, vari- 



ations in nx and T / S are highly anti-correlated, and the 
error bars on T/S do not properly reflect the ability of 
the experiments to detect the tensor signal. Other choices 
of parameters, e.g. the ratio of tensor power to scalar 
power at £ = 50, should yield smaller errors. Appendix |A| 
discusses this situation. One can write the general com- 
bination of riT and T/S as A = T/S + xut- One then 
varies x to minimize the errors on X. At this minimum, 
the uncertainties on X are uncorrelated with those on ut; 
moreover, the marginalized errors on X (with ut vary- 
ing) are identical to those on T/S in the case where nx is 
held fixed. The errors on the best choice of X are shown 
in Table || and better reflect the constraint on the tensor 
signal. 

Detection Threshold — For this set of ACDM models, 
the tensor signal could be detectable with Planck at high 
significance for T/S as low as ~ 0.1 (see Tab. |J). Only 
strong tensor signals (T/S > 1) can be isolated by MAP. 
Redshift-survey information can help, especially in case 
the tensor polarization signal is obscured by foregrounds. 
Achieving the best performance depends on b eing able to 
separate th e E and B channels of polari zation (Zaldarriaga 
et al.Q997; Kamionkowski & Kosowsky 1998); for example, 
if the i?-channel is completely ignored, then the X errors 
on the T/S = 0.1 model in Table | increase by 50%. How- 
ever, much of the leverage comes from temperature data; 
if the uncertainty from other large-angle contributors can 
be removed, the tensor plateau at I « 50 can be detected. 
For example, as the optical depth in the fiducial model in- 
creases, the detection threshold for T/S drops because r 
can be better const rained from the temperature data itself 
(Zaldarriaga et al. 1997). 

Consistency Relation — If the tensor-to-scalar ratio ap- 
proaches uni ty, the tensor tilt can be measured as well 
(Knox 1995; Zaldarriaga et al. 1997). This could allow 
a test of the inflationary prediction T/S ~ — 7?ir- In 
attempting to test specific inflationary models, e.g. the 



T j S ~ 7(1— ns) relation of power-law inflation, one should 
estimate T/S and ns at the same scale. If a ^ 0, the 
scalar tilt at scales somewhat larger than /cgd will be more 
uncertain that ns(fcfid), leading to weaker constraints on 
inflation. 

6. ASSUMPTIONS 

6.1. CMB Foregrounds 

We have assumed that a number of the frequency chan- 
nels measured by the C MB experiment s will be used for 
foreground remov al (e.g., | Brandt et aL 19 94; Tcgmark & 



Efstathioul 1996; |Bcrsanclii et all 1996; ITegmarkl 1998; 



Hobson et al. 1998), leaving only the subset listed in Table 



II available for cosmology. Assuming that foregrounds can 
De eliminated to this level may be optimistic, especially 
for the polarization at the largest angular scales and for 
small angular scales in general. For the former, the cos- 
mic signal is small (see Fig. [l] inset) and must be detected 
again st potentially larger Galactic foregrounds (Keating 



et al 



1997). Large-angle temperature signals are also use- 
ful for distinguishing curvature effects but will in any case 
be severely limited by cosmic variance. On the smallest 
angular scales, point source subtraction may be insuffi- 
cient especially where the cosmic si gnal is falling d ue to 



the finite duratio n of l ast scattering (Toffolatti et aO 1998; 
Guidcrdoni ct a l. 1998; Fcgmark & de Olivcira-Costal 1998; 
Kcfregier ct"aT| l 998 ) • 



If we eliminate even more frequency channels from the 
Fisher matrix, the errors on cosmological parameters for 
the model in Table || increase only slightly. Restricting 
MAP to use only the 90 GHz channel for temperature 
data increases errors bars by less than 10%. Using only 
the 90 GHz channel for both the temperature and the po- 
larization data is a more serious loss: 50% on most pa- 
rameters and 150% on r. Essentially, the 90 GHz chan- 
nel is nearly sample-variance limited for the temperature 
anisotropies, so little is added by the lower frequencies, 
while the channels are all roughly equally important for 
the noise-dominated large-angle polarization signal. 

For Planck, retaining only the 143 GHz channel for cos- 
mology increases the error bars by ~ 10% for tempera- 
ture only and by ~ 15% for temperature and polarization 
(but 25% on T/S). Retaining only the 217 GHz channel 
does a little better (~ 7%) on temperature and somewhat 
worse (~25%) on polarization (and a factor of 3 on T/S). 
Planck's expected polarization performance doesn't quite 
reach the sample- variance limit, but a single channel of 
temperature data saturates the limit down to a beam scale 
that varies only slightly between the prime channels. 

Conversely, if one adds the channels reserved for fore- 
ground subtraction back into the cosmological Fisher ma- 
trix, the cosmological error bars don't improve signifi- 
cantly. Even if the remaining eight channels on Planck 
were assumed to give cosmological signal, including them 
would reduce the error bars in Table || by only ~ 1%! MAP 
temperature results are similarly insensitive to the 20 GHz 
and 30 GHz channels, although the polarization signal can 
be helped (40% improvement on t, less on other quanti- 
ties). Reserving these channels for foreground subtractions 
thus comes at little cost for cosmology. 

6.2. Gravitational Lensing 



11 



Table 7 



Marginalized Errors as Function of k, 









Linear P(k) 




Smooth P(k) 


Quantity 


MAP 


0.05 


0.1 


0.2 


0.4 


0.1 


0.2 


0.4 


h 


0.22 


0.096 


0.029 


0.012 


0.009 


0.100 


0.089 


0.085 




0.24 


0.098 


0.036 


0.016 


0.014 


0.10 


0.100 


0.099 


h 

1 Lffi 1 1 


U.U ( o 


U.UvSo 


0.018 


0.011 


u.uuy 


0.036 


0.036 


0.036 


n A 


0.19 


0.081 


0.042 


0.024 


0.021 


0.088 


0.087 


0.086 


n' K 


0.055 


0.030 


0.015 


0.010 


0.009 


0.029 


0.022 


0.020 


ln(Umh 2 ) 


0.095 


0.094 


0.077 


0.054 


0.049 


0.088 


0.071 


0.065 




u.uuu 




0.050 


0.038 


u.uou 


0.058 


0.049 


0.044 


m v (eV) oc f!„h 2 


0.58 


0.55 


0.33 


0.17 


0.17 


0.43 


0.43 


0.42 




0.048 


0.048 


0.040 


0.029 


0.027 


0.045 


0.039 


0.037 


a 


0.018 


0.018 


0.015 


0.012 


0.007 


0.018 


0.013 


0.008 


lnP s (fc fld ) = In A| 


0.43 


0.43 


0.36 


0.26 


0.24 


0.41 


0.35 


0.33 


T/S 


0.18 


0.18 


0.16 


0.14 


0.14 


0.17 


0.17 


0.17 


T 


0.022 


0.022 


0.021 


0.021 


0.021 


0.022 


0.021 


0.021 


lncrg 


0.14 


0.12 


0.070 


0.038 


0.034 


0.12 


0.11 


0.102 


ln(o- 50 /o- 8 ) 


0.15 


0.066 


0.028 


0.011 


0.010 


0.030 


0.012 


0.011 


ln/3 




0.090 


0.068 


0.053 


0.045 


0.094 


0.093 


0.091 


In 6 




0.22 


0.087 


0.041 


0.035 


0.24 


0.22 


0.21 



NOTES.— All models have Q m = 0.35, fl B = 0.05, H A = 0.65, h = 0.65, and T/S = 0. n T = and cannot vary. All errors are 1 - <r. CMB 
data for MAP with temperature and polarization information included on all columns. Linear P(k) means using actual linear power spectrum. 
Smooth P(k) columns use a fitting formula that captures the break at the sound horizon but eliminates the baryon oscillations in P(k). 



In this paper, we ignore the effects of gravitational lens- 
ing on the CMB power spectrum. Gravitational lensing 
from large sca le structure smooths out the hi gh multipoles 



Cole & Efs- 



of the CMB r|Blanchard fc Schneiderj 1987, L 
tathic|u| 1989; Seljak 1996a) in a way that is dependent on 
the amplitude of the matter power spectrum at low red- 
shift. Adjusting fl\ and Qk to keep a constant angular 
diameter distance allows large changes in this amplitude 
there by breaking the angular diameter distance degener- 
acy ( Mctcalf fe Silk| 1997). However, the effect is small 
and even Pl anck cannot use it to at tain good estimates on 
Vt m and ft ( Btompor fc Efstathiou 1998). We neglect it 
here since the degeneracy is broken much more effectively 
by additional information from redshift surveys or other 
sources. 



0. 6. Varying k n 



We have assumed that the galaxy power spectrum fol- 
lows linear theory on scales longward of a wavenumber 
t mal . Non- linear evolution on smaller scales may obscure 
the cosmological information in the linear power spectrum. 
To be conservative, we have neglected all cosmological in- 
formation on smaller scales. 

Baryonic features in the linear power spectrum are of 
particular importance for cosmological parameter estima- 
tion, yet they are washed out at second -order in pcr turba- 



EHT found 



tion theory (e.g. Jain & Bertschinger 1994). 
that fluctuation levels up to fc J P(k) /2tt 2 « 0.5 preserved 
the features for the model of Table 0. Simulations give sim- 



ilar results ( Meiksin et al. 1998). Hence, for this model, 
with either CMB or cluster abundance normalization, one 
may expect linear theory to apply to fc max = 0.1ft Mpc -1 
but not much beyond. This is enough to see the first acous- 
tic oscillation. 

Because the appropriate value of fc max depends on nor- 



malization, we show the results for our fiducial ACDM 
model as a function of fc max in Table 0. As expected, 
the results do depend strongly on fc max . In particular, 
as fc max increases from 0.05ft Mpc -1 to 0.2/iMpc , the 
errors on ft and related quantities drop sharply. In this 
model, the break in the power spectrum occurs at about 
0.03ft Mpc" 1 and the first peak is at 0.08ft Mpc -1 (Fig. |). 
Little information is added between fc max = 0.2ft Mpc -1 
and fc max = 0.4ft Mpc -1 here. On these small scales, the 
acoustic oscillations have been damped away even in linear 
theory. 

While our choice of fc max may be appropriate for bary- 
onic features, other aspects of cosmology may not be so 
fragile. For example, the broadband level of the linear 
power spectrum, affected by ns and a, could potentially 
be reco nstructed from the quasi-linear regime ( [Peacock & 



Dodds 1994). As a means of exploring this, we consider an 
alteration to our treatment of the matter power spectrum, 
replacing the true linear power spectru m by a smoothed 
spectrum based on the fitting formula of Eisenstein & Hr 
(1998b). This formula preserves the sharp break at the 
sound horizon (fc w 0.03ft Mpc -1 here) but ignores all 
the smaller-scale oscillations. The resulting marginalized 
errors are considerably worse, typically equivalent to us- 
ing the true power spectrum out to fc max « 0.05ft Mpc -1 . 
Moreover, the results show little improvement as fc max in- 
creases from 0.1ft Mpc" 1 to 0.4ft Mpc -1 . In other words, 
within the constraints available through CMB satellites, 
little additional cosmological information is gained by de- 
tection of a featureless matter power spectrum. Only de- 
tection of the acoustic oscillations or spectral breaks (e.g. 
massive neutrinos) produces significant improvements. An 
exception to this conclusion is the error bar on a, the run- 
ning of the scalar tilt, which continues to improve as fc max 
increases. 
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Table 8 



Marginalized Errors as Function of Parameter Space. 



Quantity MAP (TP) MAP + SDSS 



h 


0.22 


0.22 


0.064 


0.052 


0.063 


0.048 


0.022 


0.029 


0.029 


0.024 


0.014 


0.024 


0.014 


0.012 




0.24 


0.23 


0.098 


0.081 


0.097 


0.074 


0.034 


0.036 


0.034 


0.036 


0.020 


0.036 


0.019 


0.016 


\hjflh 


0.078 


0.076 


0.042 


0.035 


0.041 


0.031 


0.015 


0.018 


0.015 


0.015 


0.008 


0.015 


0.008 


0.007 


o. 


0.19 


0.18 


0.098 


0.081 


0.097 


0.074 


0.034 


0.042 


0.035 


0.036 


0.020 


0.036 


0.019 


0.016 


n K 


0.055 


0.054 












0.015 


0.012 












ln(Q m /i 2 ) 


0.095 


0.069 


0.086 


0.073 


0.084 


0.065 


0.035 


0.077 


0.059 


0.033 


0.020 


0.033 


0.018 


0.018 


\n(n B h 2 ) 


0.060 


0.042 


0.055 


0.054 


0.053 


0.050 


0.025 


0.050 


0.035 


0.028 


0.028 


0.028 


0.028 


0.023 


m u (eV) <x f2„/i 2 


0.58 


0.53 


0.58 




0.56 






0.33 


0.30 


0.31 




0.31 






Y P 


0.020 


0.020 


0.020 




0.020 






0.020 


0.020 


0.020 




0.020 






ns(fcfid) 


0.048 


0.028 


0.045 


0.039 


0.045 


0.037 


0.014 


0.040 


0.024 


0.021 


0.018 


0.021 


0.018 


0.012 


n S (Ho) 


0.17 


0.15 


0.16 


0.14 


0.045 


0.037 


0.014 


0.14 


0.12 


0.14 


0.14 


0.021 


0.018 


0.012 


a 


0.018 


0.018 


0.016 


0.016 








0.015 


0.014 


0.014 


0.014 








lnP$ (fe fld ) = lnA| 


0.43 


0.25 


0.40 


0.35 


0.40 


0.34 


0.12 


0.36 


0.21 


0.18 


0.15 


0.18 


0.15 


0.101 


lnP^ffo) 


0.71 


0.29 


0.71 


0.57 


0.61 


0.51 


0.18 


0.61 


0.23 


0.49 


0.46 


0.27 


0.24 


0.15 


T/S 


0.18 




0.17 


0.16 


0.16 


0.15 




0.16 




0.12 


0.12 


0.087 


0.085 




T 


0.022 


0.021 


0.022 


0.021 


0.021 


0.020 


0.020 


0.021 


0.020 


0.021 


0.021 


0.020 


0.020 


0.020 


lncr 8 


0.14 


0.14 


0.13 


0.056 


0.12 


0.052 


0.035 


0.070 


0.069 


0.062 


0.027 


0.062 


0.027 


0.027 


ln(ff5o/o-8) 


0.15 


0.15 


0.074 


0.059 


0.069 


0.057 


0.033 


0.028 


0.028 


0.027 


0.017 


0.027 


0.015 


0.015 


ln/3 
















0.068 


0.058 


0.047 


0.046 


0.045 


0.043 


0.038 


In ft 
















0.087 


0.087 


0.082 


0.031 


0.080 


0.031 


0.029 



NOTES.— Q m = 0.35, tt B = 0.05, fl A = 0.65, h = 0.65, and T/S = 0. n T = and cannot vary. All errors are 1 - a. CMB data for MAP 
with temperature and polarization information included on all columns. SDSS columns use CMB data and fc max = O.l/i Mpc -1 . Blank entries 
indicate that the parameter has been held fixed. 



6.4. Smaller Parameter Spaces 

We next consider the effects of removing parameters 
from our model space. Removing parameters in the face 
of degeneracies is often motivated by Occam's razor or on 
the grounds that external information will eliminate some 
of the options available. In addition, doing so allows us 
to explore the extent of the degeneracies and facilitates 
comparison to previous works. 

Removing parameters does not require recalculation of 
models; rather, it means that we hold their value fixed at 
the fiducial value (by a prior) while the errors on other 
quantities are computed. In Table |^, we show several dif- 
ferent cases, where parameters held fixed are denoted by 
blank entries, both with and without redshift survey in- 
formation. The error bars necessarily decrease as degrees 
of freedom are removed. It is important to remember that 
improvements found by removing a parameter depend on 
which the remaining parameters are; in other words, state- 
ments such as "removing a is negligible" apply only to the 
particular context of our 12-dimcnsional parameter space. 

Tensors — Eliminating tensors from the model (i.e. as- 
suming T/S — 0) allows ns, fl m h 2 , and flph 2 to be better 
determined by MAP. Planck's longer lever arm allows it 
to constrain ns regardless of tensors. 

Curvature — Fixing the curvature of course breaks the 
angular diameter distance degeneracy by assumption. There- 
fore, CMB data alone can turn the location of the acoustic 
peaks into a measure of h, fl m , and other quantities. It is 
interesting to note, however, that MAP doesn't get enough 
accuracy on, e.g., fl m h 2 to keep uncertainties in the sound 
horizon from propagating into the measure of the angular 
diameter distance. Planck without polarization does only 
slightly better in this regard. 



If SDSS information is included, assuming fix = 
makes very little difference to errors on h and fl m , since 
the angular distance degeneracy is already broken. It does 
however affect other parameters that determine the peak 
locations such as fish 2 , fl m h 2 and ns(fcfid)- To see why 
this is so, consider the general case where fix is allowed 
to vary. One has two uncertain quantities, fl\ and fix, to 
map a fairly well-constrained quantity, the sound horizon, 
to the location of the peaks in Ci and P(k). This yields 
good constraints on fl\ and fix- Now, with fix — as- 
sumed, these two observations can both constrain fl\ and 
reduce the remaining uncertainties on the sound horizon. 
Hence, we see that assuming fix = reduces error bars 
on fl m h 2 , fish 2 , and ns(fcfid). 

Neutrinos and Helium — Removing fl v h 2 and Yp from 
the fix — case makes further small improvements. Yp 
affects the ionization history and free electron density so 
as to change the sound horizon and damping length, fl v h 2 
is partially degenerate with fl m h 2 . 

Running of the Tilt — Removing a makes very little dif- 
ference except on quantities that depend on extrapolating 
the initial power spectrum beyo nd the well-obser ved range, 



e.g., ns(H ) and A 2 S (see also Copeland ct al. 1998). As 



discussed in §0, with a free, the spectral tilt becomes scale- 
dependent, but there exists a scale at which the error on 
the tilt is unchanged by the removal of a. From the fact 
that the errors in Table [sj on ns(fcfid) change very little as 
a is removed, one can infer that fc n d = 0.025 Mpc -1 is close 
to this "pivot point" . With SDSS information, removing 
a helps to better determine T / S. 

Combined — For MAP alone, removing tensors on top of 
removing fix, fl v h 2 , Yp, and a makes a significant dif- 
ference (compare the last two columns of the MAP-only 
section of Table pi) . Apparently, combinations of parame- 
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ters were conspiring to hide their effects at I < 100. 

In short, stripping the parameter space down to a min- 
imum of baryons, CDM, cosmological constant, and tilt 
(with t and A 2 S controlling the normalization) makes a 
factor of 3 difference in the error bars from MAP even 
for parameters that are not associated with strong degen- 
eracies. Including SDSS reduces this dependence; in fact, 
MAP plus SDSS in the general parameter space with only 
curvature fixed performs as well as MAP alone in this mini- 
mal parameter space! Degenerate parameters are of course 
affected much more strongly. 

7. CONSISTENCY 

Even if the CMB does end up spinning a seamless tale 
of structure formation and cosmological parameters, it will 
not spell the end of cosmology. Only with stringent consis- 
tency checks from other types of cosmological data can we 
be confident that we have eliminated systematic effects in 
the data and its analysis. Moreover, parameter estimation 
is only as good as its underlying parameter space; testing 
CMB conclusions against other cosmological probes is an 
important way to search for unrecognized physical effects. 

One of the most important data sets for this task is the 
galaxy power spectrum from redshift surveys. While we 
have focused in this paper on the ways in which such sur- 
veys can complement CMB data to improve parameter es- 
timation, this has assumed that the measured power spec- 
trum is consistent with the locus of allowed models from 
the CMB. In fact, with the precision of upcoming surveys, 
this is not guaranteed: there are many possible spectra 
that will simply be inconsistent with our understanding 
of cosmology from the CMB. Attributing the discrepancy 
to non-linear clustering or galaxy bias has consequences 
that are testable within the survey data; it is not clear 
that the discrepancy will be resolved in favor of the CMB. 
Explanation s involving alterat ions t o the dark matte r sec- 
tor (see e.g. [Turner fc White] 1997; |Caldwcll ct~aT| 1998; 
Hu 1998 for recent suggestions) could severely modify the 
implications of the CMB for both low-redshift cosmology 
and particle physics. 

It is also possible that CMB data and the matter power 
spectrum will tell a consistent story, while the true nature 
of the universe is subtlely otherwise, causing other cos- 
mological measurements to differ from predictions. Exotic 
late-time equations of state for smooth components are an 
example of this situation. Therefore, in the following, we 
consider consistency checks with other cosmological data 
sets. 

7.1. Hq, fl m , and Acceleration 

A rich area for consistency checks is the sector of clas- 
sical cosmology: Hq, 0, m , Q\ (and implicitly Qk = 
1 — O m — 51a)- As explained above, CMB anisotropics 
suffer from a severe degeneracy here, but a large variety of 
other precision measurements are available to clarify the 
ambiguity and provide consistency checks. There are a 
number o f candidates: the m atter power spect rum, super- 
novae la (Perlmutter et al. 1998; Riess et al. 199 8), and 
direct measure s of Hp [ejgT Freedman ct al] 1998; Bland 



ford fc Kundic 1996; JCooray_et al. 1998) seem particularly 
promising. Figures W and |5| shows how the combination 
of each of these with CMB data yields a small region in 
parameter space. The overlap of these three regions would 
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Fig. 4. — Constraint regions in the Q m -h plane from various com- 
binations of data sets. MAP data with polarization yields the ellipse 
from upper left to lower right; assuming the universe flat yields a 
small region (short-dashed line). SDSS (fcmax — 0.1/iMpc -1 ) Rives 
the vertical shaded region; combined with MAP gives ills small filled 
ellipse. A projection of future supernovae la results (Tegmark et al. 
1998a, middle prediction) gives the solid vertical lines as bounds; 
combining this with MAP yields the solid ellipse. A direct 10% 
measurement of Ho gives the long-dashed lines and ellipse. All re- 
gions are 68% confidence. The fiducial model is the f! m = 0.35 
ACDM model. 
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Fig. 5. — As Figure ^, but for constraints in the Cl m -^A pl aI ¥r 
Lines and shadings are unchanged in meaning. Unlike Figure U, 
assuming the universe flat (short-dashed line) yields a line, not an 
ellipse. SDSS-only constraints are not shown. 



be a highly non-trivial test of cosmology. Direct measure- 
ments of f2 m , for example from M/L (e.g. Garlberg et al 



1997b) or clus ter evolution (e.g. Carlberg et al. 1997a; 



Bahcall et al 



1997), would also be powerful and con- 
straints on Q\ from gravitational lensing (e.g. Kochanek 
1996) provide a consistency test. 

If the smooth component of missing energy is more com- 
plicated than a cosmological constant or curvature, then 
the constraints from CMB plus the galaxy power spec- 
trum can still be compared to results from classical Hq 
programs. However, the supernovae then yield a mea sure- 
mcnt of this exotic equation of state ( Garnavich ct al. 1998; 
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Fig. 6. — Allowed region in the b-r plane for the ACDM 
model. The 68% confidence region is shown for MAP+SDSS, with 
and without polarization, and Planck+SDSS, with polarization. 
U = O.lftMpc" 1 . 



Hu et aL 1998) 



7.2. Normalization, Reionization, and Bias 

Another fruitful area for consistency checks involve pa- 
rameters associated with the amplitude of fluctuations in 
the CMB and galaxy power spectrum. 

Several sources of ambiguity in the interpretation of 
the relative amplitudes today should be resolved with the 
larger dynamic range of upcoming CMB and galaxy data 
sets. Tensors affect the COBE normalization and the spec- 
tral tilt adjusts its extrapolation to smaller scales, but the 
new satellites will focus on smaller angular scales for a di- 
rect comparison to the scales probed in redshift surveys. 
Massive neutrinos affect the matter power spectrum far 
greater than the CMB and hence change the relative nor- 
malization between these two. Fort unat ely, they should 
also have detectable signatures (see 



5.4) 



Within our parameter space, there are three remaining 
factors that alter the observed relative amplitude of the 
CMB anisotropies and the galaxy power spectrum. First, 
the growth factor between recombination and the present 
day depends upon Q\ and £Ik- Second, galaxy bias al- 
ters the normalization of the galaxy power spectrum rela- 
tive to that of the matter (and CMB). Third, reionization 
suppresses CMB anisotropies for a given level of potential 
fluctuations. 

Because the addition of redshift survey information 
breaks the angular distance degeneracy, the growth fac- 
tor will be accurately known. As discussed in the previous 
section, there are many consistency checks to ensure the 
validity of this statement. 

This leaves the bias and reionization optical depth 
strongly covariant, as shown in Figure |(| Without po- 
larization information to determine t and other dynam- 



ical measurements to determine b, neither will be well- 
determined. 

With a polarization detection of the large-angle reioniza- 
tion signal, however, the error bars on r plummet, yielding 
b and er 8 to ~ 5% fractionally. Because polarization mea- 
surements arc subject to many systematic effects, consis- 
tency checks here are particularly important. With Q TO 
measured precisely, any means of constraining (3 or as will 
provide a consistency check on t his measurement, e.g. red- 
shift distortions (see Hamilton 1997 for a review), pecu- 
liar velocity catalogs ( Strauss fc Willick 1995 and refer- 
ences within) , an d cluster abundances (White et al. 1993b; 



Viana & Liddle 1996 



Eke et al. 



1997; penj 1997Y Addi 



tionally, if these measurements can produce limits on (3 
or below the values given in the Tables (<; 5%), then this 
extra information can place leverage on other parameters 
such as h and O m /i 2 once r has been fixed by the polar- 
ization detection. 

Conversely, if large-angle CMB foregrounds prevent us 
from extracting the reionization signal, precision measures 
of b or (is would sharply constrain r, thereby requiring 
our cosmological models to ionize the universe at some 
particular redshift. If this redshift is uncomfortably high 
or low, it may challenge our cosmological assumptions, for 
example, our extrapolation of P$(/c) to yet smaller scales. 

8. CONCLUSION 

We have presented Fisher matrix calculations of the 
cosmological information obtainable with upcoming CMB 
satellite missions and large redshift surveys. We have used 
a considerably larger parameterization of adiabatic CDM 
than previous works involving polarization data or redshift 
surveys. Within this space, we have conducted several pa- 
rameter studies, including variations in fls, and T/S. 
In a number of cases, we find softer degeneracies and hence 
larger errors bars than prior work. We attribute these 
discrepancies to artificially broken degeneracies caused by 
numerical subtleties such as use of one-sided derivatives 
and unfortunate step sizes. 

The primary purpose of this paper has been to explore 
the ways in which galaxy power spectrum data can provide 
complementary information to the CMB, thereby signifi- 
cantly reducing error bars. In the language of degenera- 
cies, we seek places in which large-scale structure offers the 
means to break CMB degeneracies and vice versa. The 
most important example of this is the angular-diameter 
distance degeneracy of the CMB, which causes uncertainty 
on h, ft m , £1a, and even as- The presence of baryonic os- 
cillations in matter power spectrum offers a robust way 
to break this degeneracy. Redshift-survey data also helps 
determine the mass of the heaviest neutrino, especially if 
it exceeds 1 eV. Although the projections of the Planck 
satellite with polarization suffice to do most everything 
else, less sensitive CMB experiments can benefit from red- 
shift surveys even on measuring quantities such as fl m h 2 
and fi^/i 2 . 

While CMB and galaxy power spectra themselves offer 
a myriad of consistency tests, they also provide a baseline 
for tests against other cosmological measurements. With 
the detection of the large-angle polarization signal due 
to reionization, the two data sets together yield a ~ 5% 
measure of linear galaxy bias, to be compared to that ob- 
tained from peculiar- velocity and redshift-distortion meth- 
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ods. Similarly, supernovae distance measurements and 
Hubble constant measurements provide a fertile set of 
cross-checks on the breaking of the angular-diameter dis- 
tance degeneracy. 

It is important to remember that the error estimates pre- 
sented here reflect the statistical leverage available within 
these data sets. We have not attempted to produce a data- 
analysis pipeline that would address the horde of obstacles 
that stand between the raw data and these cosmological 
inferences. Considerable effort has been channeled toward 
the development o f su ch methods. Foreground removal 
was mentioned in 



3.1 



Pipeline methods have b een de- 
velop ed and imple mented fo r CMB mapmaking (Wright 
1996; |Wright| 1996; |Tcgmark| 1997), CMB power spectrum 
extraction "( [Tegmark| 1997c; |Bond et alj 199 8; |Oh et al] 
19 98) and galaxy power spectrum estimation (Tegmark et 
al. | 19f )8b). All of these are lossless in the sense that they 
retain all the cosmological information quantified by the 
Fisher information matrix, but it is likely that many as- 
pects of the problem will only be revealed once the data is 
in hand. Approaching the performance described by anal- 
yses such as those in this paper will be one of the primary 
goals (and motivations for) of all this hard work! 

In summary, while CMB satellite missions can provide 
marvelous cosmological information, they do not make 
other methods obsolete. Additional precision measure- 
ments will be needed both to break the parameter degen- 
eracies of the CMB and to test for consistency in the face 
of systematic errors and additional physical effects. Only 
with these precise comparisons can one build a secure cos- 
mological model. 
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A.l. Proof 

We assume a n-dimensional parameter space with in- 
dependent variables p\, . . ., p n . Let the primary quan- 
tity be pi and the m — 1 quantities (m < n) that we 
want to combine be p%, .. . ,p m - We seek a new quantity 
X = f(pi,p2, . . . ,Pm) that has df /dpi — 1 and minimum 
errors. The former requirement is to prevent simple rescal- 
ings of the variable p%] X and p\ will in this sense have 
the same scale. 

In the Fisher matrix formalism, all that matters is the 
gradient of X with respect to the independent variables. 
We write Wj = df/dpj and arrange them as a vector w; 
by construction, wi — 1 and Wj = for m < j < n. Then 
the smallest attainable variance of X will be 



(Tv = w Cw, 



(Al) 



where C is the inverse of the Fisher matrix F. We seek 
the set of W2, ■ ■ ., w m that minimize a\. Clearly only the 
submatrix C involving the first m rows and columns of C 
can be involved. Using Lagrange multipliers, one readily 
shows that the minimum is achieved at 



(C- 



(C- 




r 1 3 < ™, 
in 

m < j < n. 



(A2) 



The minimum value of a\ is then (C -1 )^ 1 . 

If one replaces the independent variable p\ by A, then 
the transformed matrix C will have C\j = for j = 
2, . . . , m. In other words, X is uncorrelated with the vari- 
ables P2t ■ ■ iPm- This means that X will have the same 
errors regardless of whether p 2 , ■ ■ ■ , p m vary or are held 
fixed. If they are held fixed, then X and p\ are identical. 
Hence, the error on X with P2, ■ ■ ■ ,p m varying is the same 
as the error of p\ with them held fixed. 

Note that in the case of m = n, C* -1 = F. Hence, w 
is just the renormalized first column of the Fisher matrix, 
and the minimum variance of X is simply 1/Fn. Aficiona- 
dos will recognize this as the variance of p\ if all other 
variables are held fixed. 

In the case of m = 2, where we wish to combine p\ 
with one other variable pi, we obtain the special case that 
X = pi -P2C12/C22 and a\ = C u -C^Cw X and p 2 
are uncorrelated. 



APPENDIX A 

MINIMIZING ERRORS BY REDUCING COVARIANCE 

At several points in the paper, we are concerned with 
how the uncertainties on a particular quantity depend 
upon the uncertainties on other related quantities with 
which it is correlated. Here we show that a particular lin- 
ear combination of the correlated quantities has the min- 
imum variance and is uncorrelated with the other param- 
eters. This variance is furthermore equal to that of the 
original quantity under the assumption that the correlated 
parameters have been fixed by priors. We examine opti- 
mal combinations of the scalar tilt ns and running of the 
tilt a as examples. The case of the tensor-scalar ratio T/S 
and tensor tilt nx was discussed in § [3. 5] , 

7 http://arct urus.mit.edu/~matiasz/CMBFAST/cmbfast. html 



A. 2. A tilt example 

An interesting example of this concerns ns and a. We 
define our independent variable to be ns(fcfid) and a. For 
most values of fcfid, these two will be correlated. What is 
special in this case is that the minimum-error combination 
A is a physically-motivated quantity, namely the tilt ns 
at some new scale fc p i vo t = fcfid exp(— C nsa /C aa ) (see eq. 
Jill). As shown above, ns(k p i vot ) is uncorrelated with a 
and has the same error with a varying as ns on any scale 
would have if a were held fixed. 

The pivot scale fc p i vo t at which the error on ns(fc p i vo t) 
is minimized depends on the experiment and on the fidu- 
cial model. For our ACDM model and SDSS alone, fc p i vo t 
is 0.024 Mpc" 1 (0.088 Mpc" 1 ) for fc max = 0.1/iMpc" 1 
(0.2/i Mpc" 1 ). For CMB data (with SDSS to 0.1/iMpc^ 1 
in parentheses), the scales in Mpc -1 are 0.034 (0.036) for 
MAP(T), 0.018 (0.020) for MAP(TP), 0.084 (0.070) for 
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Table 9 



List of Independent Parameters. 



Notes 



Quantity 


Step size 


tlmh 2 


±5% 


U B h 2 


±5% 


n v h 2 


±10% 


n A 


±0.05H m 


n K 


A^ m = -o.in 


T 


±0.02 


Y v 


±0.02 


"s(fcfld) 


±0.005 


a 


±0.005 


T/S 


Exact 


iij- 


±0.01 


4 2 


Exact 


P 


Exact 



Equivalent to m„ 



See S B.2.1 



Prior of ±0-02 
See eq. (IA) 
See eq. 

Only if T/S + 0. 



NOTES. — Step sizes are for our ACDM model; those listed as per- 
centages are fractions of the fiducial value. 



Planck(T), and 0.029 (0.027) for Planck(TP). As expected, 
Planck adds more small-scale sensitivity, while polariza- 
tion adds more large-scale sensitivity by resolving the low- 
l degeneracies. Our choice of fcgd = 0.025 Mpc -1 is close 
to the desired spot, and the errors grow only to second 
order in ln(fc) away from the minimum: 



is(fcpivot) 



ln 2 (fc/fc pivo t)CT^. 



(A3) 



APPENDIX B 
NUMERICAL METHODS 

B.l. Derivative Methodology 

As described in §|2|, the calculation of the Fisher matrix 
reduces to manipulation of derivatives of the various power 
spectra with respect to cosmological parameters. How- 
ever, the near cancellation of certain linear combinations 
of derivatives leaves the Fisher matrix nearly singular. As 
the larger error bars are themselves inverses of the smaller 
eigenvalues, it is critical to prevent numerical effects from 
perturbing these small eigenvalues. Hence, constructing 
and manipulating the derivatives requires care, lest a pa- 
rameter degeneracy be broken by numerical effects and 
yield an overestimate of the experiment's ability to mea- 
sure the associated parameters. We find that our treat- 
ment of certain derivatives reveals significantly softer di- 
rections in parameter space than found by previous works 
(see Appendix |c]). Therefore, we will describe our meth- 
ods in some detail. 

B.l.l. Two-sided derivatives 

Wherever possible, we take two-sided derivatives. Writ- 
ing f(p) for the dependence of either lnC^ or lnP(fc) on 
some parameter p with all others fixed, this means that 
we approximate f'(p) by 

/'(,) . gg±Mzi(£z^) (B1) 



for some small step size Ap. This is exact to 2 nd order 
in Ap, whereas approximating f'(p) with the one-sided 
difference [f(p + Ap) - /(p)]/Ap « f'(p + Ap/2) is only 
good to first order; moreover, the latter corresponds to 
an accurate estimate of the derivative at a slightly shifted 
parameter value p+Ap/2. This is critical when perturbing 
parameters that change the locations of the Doppler peaks, 
because the various derivatives will no longer be in phase! 
If the step size is sufficiently large, these phase shifts will 
break the parameter degeneracies. 

One also should avoid differencing models that are cal- 
culated with different numerical techniques, as these can 
cause discontinuous results as one adjusts the independent 
variable. With CMBfast v2.3.2, this situation occurs for 
geometrically flat versus open models and for models with 
differing numbers of species of massive neutrinos. 

Because of the desire to use two-sided derivatives, we 
take non-zero fiducial values for non-negative quantities 
such as r = 0.05 and f2„/f2 m = 0.05 unless otherwise 
noted. Perturbing around fi„ ^ also allows us to use 
one species of massive neutrinos in all cases. 

B.l. 2. Derivative step sizes 

What is the best choice of the step size Ap for con- 
structing the derivative in equation (Bl)? We use CMB- 
fast v2.3 .2 (inc luding the bug fixes o f v2.4.1; [Bcfjak fc Zal- 
darriaga 1996; |Zaldarriaga fe Seljak 1997) for our numer- 



ical derivatives and find that this version of CMBfast has 
random numerical noise at a level of a = 5Ce/Cg ~ 10~ 4 
for geometricall y flat models and 10~ 3 for open models.^ 
As described by Press et al. (1992), the optimal compro- 
mise between numerical noise and higher-order Taylor se- 
ries terms occurs for Ap ~ V 1/3 Pc, where p c indicates the 
characteristic scale on which p varies. Hence, if p c ~ p, 
we should use fractional steps of roughly 5%, yielding an 
accuracy of ~0.2%. Note that a one-sided derivative with 
a 5% step would yield a truncation-dominated accuracy of 
~5%, which is considerably worse. 

Since this estimate is quite crude, we performed a series 
of numerical experiments with different step sizes before 
arriving at the choices in listed in Table ^. We have tested 
that our answers change by less than 10% when using half 
the listed steps. Convergence is better as degeneracies are 
lifted by complementary information or by reducing the 
parameter space via priors. Only for the open fiducial 
model do we see some lack of convergence due mainly to 
the larger noise associated with open models. The step 
sizes on ns, fiT, and a could easily have been signifi- 
cantly smaller, although reducing them to 10 -4 changes 
the marginalized errors by less than 1%. 

B.2. Curvature 

Curvature is known to be strongly degenerate with other 
effects, such as a cosmological co nstant, that alte r the 
late-tim e evolution of the universe (Hu fc Sugiyama] 1995 



VI.B2; |Bond et alj 1997; |Zaldarriaga et alj 1997) Co 



binations of changes in these parameters that hold the an- 
gular diameter distance to last scattering fixed will leave 

8 This does not apply to variations in the initial power spectrum, 
as the anisotropics generated by different wavenumbers are weighted 
and combined at double precision. The quoted levels are for I < 100 
and may increase beyond that. 
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1000 



Fig. 7. — (a) Derivatives of Cti with respect to at constant 
Qa (solid line) and with respect to Qa at constant Qk (dashed 
line). The former is a one-sided derivative formed by differencing 
models with Q[f = 0.01 and Q[f = 0.003. (b) Derivative of Cti 
with respect to Six at constant d^, constructed in two different 
ways, (solid line) result from differencing two models with Qa and 
QfC altered so as to keep dj\ fixed, (dashed line) the appropriate 
linear combination of the curves in panel (a) ; the poor cancellation 
is due to the curvature derivative being one-sided and therefore out- 
of-phase. The glitch at I »s 500 in the solid curve is numerical noise 
and the reason we need to truncate this derivative. 



the acoustic peaks unchanged (given a proper choice of 
normalization). Differences in the decay of potentials at 
late-times will alter the large-angle anisotropy; this alter- 
ation of the ISW effect is the primary wa y to brea k this 
degeneracy in the absen ce of lensing effects (Hu ct al. 1997; 
Stompor fe Efstathioul 1998). 

When faced with a strong degeneracy, it is numerically 
desirable to change variables in parameter space by in- 
troducing a parameter in the degenerate direction. This 
has the advantage of containing all numerical problems in 
a single (very small) derivative, where the physical effects 
can be easily understood and distinguished from numerical 
errors making the derivative artificially large. We adopt 
this approach for the angular distance degeneracy, as de- 
tailed in the following three subsections. 



B.2.1. Differentiating with (1a fixed 

Including variations in curvature around a flat model 
presents a particular concern because no closed version of 
CMBfast is available. Hence, the derivative with respect 
to curvature is necessarily a one-sided derivative. More- 
over, the open version of the code is noisier than the flat, 
which forces one to take larger step sizes for the deriva- 
tives. Since adding curvature causes a significant shift in 
the location of the CMB peaks, large step sizes cause the 
curvature derivative to be out-of-phase with the deriva- 
tive. This is displayed in Figure ^ the dashed line in the 
bottom panel shows that the two derivatives in the top 
panel cannot be exactly cancelled. 



We work around this problem by altering our coordinate 
basis so that the only derivative in a direction with non- 
zero curvature has a convenient physical property. The 
derivative of J7a at constant fix may be done in the usual 
manner. But for the derivative of fix at constant Q,a, 
we transform to the basis of fix and angular diameter 
distance <1a- The latter is defined as 



1 



sinh 



dz 



(B2) 



e{z) = n A + Q K (1 + zf + a m (i + zf + Q R (1 + zf, 

where flu is the d ensity of radia tion and is the redshift 
of recombination ( Hu fc White 1997a). We treat massive 
ne utrinos according to the approximation of Appendix A 
of Hu & Eiscnstcin (1998). The desired curvature deriva- 
tive is then 



d 



an 



( d 



\dfl 



K 



( Ma 

\on K 



( 9 



I dd A \ 



(B3) 



In flat models, the first term is constructed by a one-sided 
finite difference in curvature corresponding to Af2 m = 
— 0.1f2 m . The second term consists of a numerical prefac- 
tor times the derivative that we have already calculated. 
The first term is special because with the proper choice of 
normalization it is a derivative in the direction of constant 
Doppler peak location, shape, and height. This means 
that aside from an ISW signal at low £ and lensing effects 
at high £, the derivative vanishes, allowing one to isolate 
the numerical noise. Also, since the peak structure of both 
derivatives enters through the same d/dfl^ at constant fix 
term, the two derivatives are guaranteed t o ca ncel up to 
the treatment of the first term of equation ( B3 ) . 



B.2.2. Differentiating with fixed high-z normalization 

To properly implement the technique of the last sec- 
tion, one must choose a normalization convention that is 
independent of Q\ and fix- This condition is satisfied by 
keeping unchanged the physical situation at high redshift. 
It is implemented by holding fixed the amplitude of the 
scalar gravitational potential at high redshift on some co- 
moving reference scale fc n0 rm in Mpc -1 that is outside the 
horizon at the redshift in question. We use the low-redshift 
growth function to shift the scalar potential back to the 
pure-matter-domination phase (ignoring radiation). This 
produces the normalization 



0) = A\ 



k 



iid 



1 + 



2 \ 2 

Kll \ tj-4 i n 2 



1-2 



gr' 



(B4) 

(ng = 1) in other words, we choose A| as an independent 
variable and therefore renormalize the calculated (COBE 



normalized) P(k) and CV by Equation (B4) to hold A s 



constant when taking all ot her deriy atives. Here D gr is 
the growth function integral ( Peebles] 1980) 



1 3 / 2 da 



(f2 m + fl K a + fi A a 3 ) 3 / 2 



(B5) 



Note that in the Harrison-Zel'dovich case, the normaliza- 
tion is independent of fc nor m m the large-scale limit; we 
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choose a normalization scale of k~^ lm — 3000 Mpc. If 
Us 7^ 1, it is important to pick fc norm to be independent 
of h. Since the high-redshift ratios of the various types of 
matter and radiation are held fixed when perturbing J7a 
and Ox, it doesn't matter whether we normalize to the 
matter-domination or radiation-domination potential. 

B.2.3. Noise clipping 

Even with the derivative technique of the last two sec- 
tions, the numerical noise of CMBfast v2.3.2 breaks the 
angular diameter distance degeneracy. Fortunately, with 
this choice of derivatives, the true signal is we ll- localized 
in I. Neglecting gravitational lensing (see §5.2), the only 
signal in the curvature derivative at constant oIa is the 
\ow-£ IS W effect. This effect scales as Q ~ f s " 4 asymp- 
totically (|Hu fe White] 1996). 

Since the goal of our analysis is to make robust and 
conservative estimates for how accurately cosmological pa- 
rameters can be measured, we assume that no useful cos- 
mological information can be extracted from lensing effects 
or from ISW contributions at large £. For £ > £ c u p = 30, 
we therefore set C' e — (£ c \ip/£) 4: ~ ns C e for this partic- 
ular derivative, thereby throwing away all of the high-£ 
numerical noise that would otherwise artificially break the 
degeneracy above some £-cutoff ^ c iip- This numerical noise 
is displayed as the solid line in Figure ^ the glitches at 
£ w 500 dominate the breaking of the degeneracy for ei- 
ther satellite. 

Hence, while our error bars on Q\ and £Ik for Planck 
alone (especially with polarization) are overestimated due 
to the neglect of lensing, we protect the softness of this de- 
generacy against the numerical problems that would oth- 
erwise dominate. 

B.2.4. Effects of il K and tt A on P{k) 

Q.K and Qa have degenerate effects on P(fc) as well, and 
so we must take care in constructing these derivatives. For- 
tunately there is an analytic solution. The matter power 
spectrum can be decomposed into an initial power spec- 
trum whose time dependence reflects only the physics on 
the largest scales times a transfer funct ion that incorpo- 
rates the effects of causal physics (e.g. Eisenstein fc Hrj 
1998ab). VLk and Vt\ shift the the z = normalization of 
the initial power spectrum relative to the level of CMB an- 
isotropics by altering the growth function D gr [eq. (B5)]. 
Meanwhile, the transfer function is independent of fix and 
f^A if it is measured in real space instead of redshift space. 
Hence, Hk and Ha enter this piece only through their ef- 
fect on h. We can rewrite this as a derivative with respect 
to k: 



rflnP _ <iln,D gr 1.2 



1 



rflnP 
din fc 



- 1 



(B6) 



This avoids having to difference two different Boltzmann 
code outputs and allows us to track the oscillations in the 
power spectrum and its d eriv ative more accurately. The 
middle term in equation (B6) comes from the derivative 
of In b, since (3 is our independent variable. The deriva- 
tive with respect to JIa has an equivalent formula. The 
assumption that the transfer function in real space is in- 
dependent of Q.k and VIa is actually violated slightly at 



small scales if massive neutrinos are imp ortant due to dif- 
ferences in the infall of the neutrinos (Hu & Eisenstein 
1998); we ignore this effect as it is tiny and primarily be- 
yond the linear regime. 

B.3. r and Normalization 

Reionization mimics a suppression of the amplitude of 
the primary CMB anisotropies on all but the largest scales. 
For small optical depths, the secondary fluctuations gener- 
ated by these late-time scatterings are small because only a 
small fraction of photons are affected and because the scat- 
tering occurs over a sufficient range of distances along the 
line of sight that small-scale perturbations are averaged 
out (Kaiser 1984). Hence, for r <C 1, the main effect is a 
suppression of Ci by exp(— 2t) for all but the lowest £. If 
the normalization of the primary fluctuations is unknown, 
then we can hide the reionization simply by increasing 
their normalization by a corresponding amount. The re- 
sulting degeneracy is difficult to break in the temperature 
anisotropy data because the rise at low I up to the original 
fluctuation level is hidden by cosmic variance. With po- 
larization data, reionization can be more easily separated 
because it produces a bump at large-angles whereas fluc- 
tuations at high redshift cann ot prod uce much large-angle 
polarization due to causality (Kaiser 1983). 

With this level of degeneracy, one should worry about 
artificially breaking the cancellation between the r and A 2 S 
derivatives. However, we find that CMBfast is sufficiently 
accurate to track this degeneracy without invoking tricks 
similar to those used for the curvature, e.g., differentiating 
with A$e~ T held fixed. 

Note that if our chosen normalization were the COBE 
normalization, we would also have to worry about degen- 
eracies between the tensor-to-scalar ratio T/S and the 
scalar normalization. By normalizing to the scalar po- 
tential fluctuations, however, we avoid this problem. 



APPENDIX C 

COMPARISON TO PREVIOUS WORK 

The results in the main part of this paper are not di- 
rectly comparable to those of previous papers due to dif- 
ferences in parameter spaces or fiducial models. In this 
Appendix, we attempt to provide as direct a comparison 
as we can for some of the results presented in these pa- 
pers. We have in all cases matched the quoted parameter 
sp aces an d experimental specifications. Except in the case 
of |Wh!tcl (1998), we have not taken care to provide an ex- 
act match to the normalization in these comparisons. All 
groups are normalizing to the COBE value for the signal 
to noise in equation (0). While normalization differences 
of a few percent are possible, these would not affect the 

parameter results beyond this small level. 

In Table |Tol we compare marginalized errors with Bone] 



et al. (1997jfor their standard CDM model. We use their 
experimental specifications for MAP and the Planck High- 
Frequency Instrument HFI (temperature only) and restrict 
ourselves to their parameter space (£Ik = a = 0). Note 
that we have by necessity used one-sided derivatives for 
r and f2„ differencing models with 0.01 and 0.001 in each 
parameter. Because the curvature is not varied, the tricks 
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Table 10 



Table 12 



Comparison to Bond et al. (1997) 



Comparison to Wang et al. (1998) 





MAP 


Planck (HFI) 


Quantity 


BET 


EHT 


BET 


EHT 


2h 


0.19 


0.40 


0.02 


0.08 


4Q A h 2 


0.49 


1.2 


0.05 


0.20 


lnQ B h 2 


0.09 


0.14 


0.006 


0.018 


4fl m h 2 


0.18 


0.38 


0.02 


0.04 




0.07 


0.09 


0.02 


0.03 


n s 


0.06 


0.10 


0.006 


0.015 


T/S 


0.38 


0.69 


0.09 


0.11 


T 


0.22 


0.34 


0.16 


0.30 


In a s 


0.28 


0.57 


0.18 


0.33 



NOTES.— Fiducial model has H 
tt A = 0, r = 0, n s = 1, T/S = 
AY P = 0.Q2 All prmr« a re 1 



quoted in 



Bond et al. 



l,h = 0.5, fl B = 0.05, Q v = 0, 
= 0, and Yp = 0.23 (with a prior of 
- a. BET columns contain the errors 
EHT columns contain the errors we 



(1997); 

find using their experimental specifications for MAP and the High- 
Frequency Instrument (HFI) of Planck. All columns are temperature 
data only. 



Table 11 







MAP(T) 


MAP(TP) 


Model 


Quantity 


WSS 


EHT 


WSS 


EHT 


SCDM 


In h 


0.052 


0.051 


0.033 


0.033 




Q A 


0.15 


0.15 


0.091 


0.097 




\nU B h 2 


0.028 


0.031 


0.020 


0.024 




In t 


2.4 


2.9 


0.39 


0.38 




ns 


0.017 


0.020 


0.0085 


0.0090 


ACDM 


ln/i 


0.066 


0.077 


0.032 


0.035 




In SI a 


0.076 


0.089 


0.037 


0.043 




lnQ B h 2 


0.044 


0.052 


0.021 


0.023 




In t 


1.3 


1.6 


0.18 


0.18 




ns 


0.035 


0.041 


0.014 


0.014 



NOTES. — Top set of numbers are for a standard CDM fiducial model 
with Q m = 1, h = 0.5, Q B = 0.05, r = 0.05, and ng = 1. Bottom 
set are for a ACDM model with Q m = 0.3, h = 0.65, O s = 0.06, 
Ha = 0.7, t = 0.1, and ns = 1- £Ik, Qv, Yp, T/S, and a are held 



fixed. All errors are 1 



WSS m. 



umns contain the errors quoted 
(1998; private communication). 



in the revised version of Wang et al. 
EHT columns contain the errors we find using their experimental 
specifications for MAP with and without polarization. 



Comparison to Zaldarriaga et al. (1997) 



source of the discrepancies 



MAP(T) 



MAP(TP) 



Quantity 


ZSS 


EHT 


ZSS 


EHT 


h 


0.092 


0.17 


0.051 


0.040 


n A 


0.53 


1.1 


0.29 


0.25 


n B h 2 


0.0010 


0.0019 


0.00061 


0.00055 


T 


0.13 


0.21 


0.021 


0.022 


n s 


0.059 


0.12 


0.031 


0.028 


T/S 


0.39 


0.77 


0.22 


0.19 


h 


0.017 


0.025 


0.016 


0.017 


n A 


0.098 


0.15 


0.093 


0.103 


fl B h 2 


0.00030 


0.00043 


0.00028 


0.00035 


T 


0.12 


0.15 


0.021 


0.021 


ns 


0.0098 


0.021 


0.0048 


0.010 



T/S 



NOTES.— Fiducial model has Q TO = 1, h = 0.5, tt B = 0.05, r = 
0.05, 715 = 1> an d T/S = 0. Qk> &.v, Yp, and a are held fixed. 
The top set of numbers is with T/S allowed to vary; the bottom set 
is with T/S fiyeH All errors are 1 — <j. ZSS columns contain the 



errors quoted in Zaldarriaga et al, (1997); EHT columns contain the 
errors we find using their experimental specifications for MAP with 
and without polarization. 



of §B.2 are not needed. We find significantly (up to a 
factor of 3) larger err ors on some parameters. 

Zaldarriaga et al.' (1997) analyzed a standard CDM 
model with a smaller parameter space. Restricting to 
this space {VIk = n w = a: = AYp = 0) and adopting 
their specifications for MAP yields the results in Table O. 
Because r ^ and massive neutrinos are excluded, no 
one-sided derivatives are required. The results agree well 
with temperature and polarization data (except for ns) 
but can be up to a factor of 2 different with temperature 
data alone. 

We s uspe ct that the numerical issues addressed in Ap- 
pendix B.l are responsible for discrepancies described 
above. Unfortunately, based on the published informa- 
tion, we cannot confirm in the above cases that this is the 



Our comparison with Wang et al. (1998) demonstrates 
the situation. We compared two fiducial models, adopting 
their parameter space (£Ik = a = ft v = T/S = AYp = 0) 
and specifications for MAP (with / s k y = 0.8). At first, 
agreement was poor, with their errors being as much as 
a factor of 5 smaller on one model. One-sided derivatives 



had been used; when Wang et al. kindly recomputed their 
results using two-sided derivatives and smaller steps on 
certain parameters, the agreement became quite good (<J 
20%), as show in Table O. We did not compare results 
for SDSS and CMB together because the treatments are 
quite different: they remove baryon oscillations from P(k) 
and apply a non-linear evolution correction. 

In another situation in which we could confirm that two- 
sided derivatives were used, we have attained even better 
agreement. We have compar ed a ACDM model with the 
preliminary results of White (1998), who uses a hierarchy 
Boltzmann code rather than CMBfast. Under a parameter 
space of ris/i 2 , fl rn h 2 , t, ns, and A 2 S , the marginalized 
error bars agree to 3% with and without polarization for 
MAP and Planck. Adding fi^ causes larger discrepancies, 
mostly along the direction of the angular diameter distance 
degeneracy; we suspect the hierarchy code is more stable 
for fl k ^ 0, but the differences should disappear once the 
degeneracy is broken by outside information. 

Degeneracies are quite important even in these smaller 
parameter s paces. As an exam ple, consider the SCDM 
model from Zaldarriaga et al. (1997). Holding r fixed 



removes the reionization-normalization degeneracy and 
leaves ri m /i 2 , B /i 2 , Q\, n s , T/S, and A 2 S to vary. This 
set of parameters in this fiducial model can be combined 
to cancel the normalization temperature derivative to bet- 
ter than 1 part in 100 over the Grange 100 — 800, mostly 
through a combination of tilt, tensors, and fl m h 2 . Hence, 
even in this small space, one needs to control the deriva- 
tives to high numerical accuracy. For example, an 0.05 
one-sided step in tilt pivoting around the Hubble distance 
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produces a 16% derivative error at I « 600 due to second- 
order temis^_These consid eration s and the comparison 
with |Wang et alj (1998) and |Whitcj (1998) lead us to con- 
clude that two-sided derivatives are crucial for achieving 
accurate answers with present-day codes. 
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